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ABSTRACT 


Work in the stabilization of monolithic Nd:YAG lasers and the 
application of these lasers to nonlinear optical frequency conversion is 
discussed. The intrinsic stability of semiconductor-diode- laser-pumped 
solid-state lasers has facilitated a number of demonstrations in external 
resonant cavity harmonic generation and stable optical parametric 
oscillation. Relative laser frequency stabilization of 0.3 Hz has been 
achieved, and absolute stability of a few hundred hertz is anticipated. The 
challenge is now to reproduce this frequency stability in the output of 
tunable nonlinear optical devices. Theoretical and experimental work 
toward this goal are continuing. 
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I Introduction 

This progress report is presented in the form of two manuscripts that have been submitted 
for publication. The first manuscript describes a presentation given in Vilnius, Lithuania in 
August of 1990. Our work in the development and frequency stabilization of 
semiconductor-diode-laser-pumped solid-state lasers was presented. Emphasis was placed 
on research in nonlinear optical frequency conversion made possible through the use of 
these lasers to generate pump radiation for the nonlinear optical devices. These lasers are 
now being stabilized to better than one hertz relative to a high-finesse resonant cavity. 
Preliminary demonstrations have shown the possibility absolute stabilization to a few 
hundred hertz using sub-Doppler spectroscopy techniques with molecular absorptions in 
iodine. 

The second manuscript describes a theoretical analysis of the tuning properties of 
monolithic doubly resonant optical parametric oscillators (DRO's). This paper will be 
published in the March issue of the Journal of the Optical Society of America B. The 
tuning properties of the DRO's are complicated, and stable operation requires that an 
number of critical tolerances be met. Monlithic DRO's that we have operated have 
displayed extremely good coherence properties, low thresholds (10 mW) for cw parametric 
oscillation, and good conversion efficiency. The potential of optical parametric oscillators 
for optical communication, remote sensing, spectroscopy, and even projection displays is 
extensive. This potential for application proves the motivation for stabilizing the DRO. 

We are currently pursuing a number of projects in these areas. Work on both 
improved relative laser frequency stabilization and absolute stabilization is continuing. The 
single-mode output-power of our injection-locked Nd:YAG laser has recently been 
upgraded to 18 W, and cw harmonic generation has be demonstrated at 6.5-W power 
levels. This is adequate power to pump even a single resonant optical parametric oscillator. 
In addition we will be studying cascaded harmonic generation to the ultraviolet and a 
number of additional DRO and SRO schemes. 
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Frequency Stable Operation of Nd:YAG Lasers 

R. C. Eckardt, T. Day, E. K. Gustafson and R. L. Byer 

Ginzton Laboratory 
Stanford University 
Stanford, California 
94305-4085 USA 

Abstract 

Semiconductor-diode-laser-pumped solid-state lasers are inherently frequency 
stable when compared with other lasers, and this stability can be improved by 
active feedback control. Monolithic nonplanar ring oscillators constructed from 
solid pieces of the laser material provide better than 10-kHz frequency stability 
over 0.1 -sec intervals. Active feedback stabilization of the cavity length has been 
used to demonstrate 0.3-Hz stabilization relative to a reference cavity resonance 
frequency. This range of frequency stability has application in nonlinear optical 
frequency conversion. 


Resonantly enhanced nonlinear frequency conversion techniques, such as extemal-resonant- 
cavity harmonic generation and optical parametric oscillation, are particularly sensitive to the 
frequency stability of the pump radiation. High stability diode-laser-pumped solid-state lasers 
permit efficient cw second-harmonic generation (SHG) and optical parametric oscillation (OPO). 
An interesting aspect of this work is an analysis of the conditions required to produce single- 
mode-pair optical parametric oscillation. The exceptional frequency stability of the diode-pumped 
NonPlanar Ring Oscillator (NPRO) compared to dye lasers, gas lasers, or semiconductor lasers 
also provides some unique advantages for active frequency stabilization. Sub-Hertz relative 
frequency stability has been achieved by locking to reference cavities, and absolute stability of a 
few hundred Hz by locking to molecular or atomic absorptions appears possible. 

Remeasurement of some nonlinear optical coefficients is a by-product that grew out of this 
work when it was found that some of the accepted values were predicting results that were 
inconsistent with experimental observations. 

1. The NonPlanar Ring Oscillator 

The solid-state monolithic NPRO 1 ’ 2 * 3 generates stable single-axial-mode output which is 
important for resonant cavity nonlinear optical frequency conversion. Semiconductor-diode-laser 
pumping of the NPRO's 4 ’ 2 provides advantages of efficiency and reduced heat loading. In 
addition, the low amplitude fluctuations of the diode-laser pump source provide improved 
frequency stability. 

A schematic illustration of a NPRO is shown in Fig. 1. Spatial hole burning in the popula- 
tion inversion is avoided by unidirectional oscillation in a ring cavity which is important for single 
frequency oscillation. The monolithic cavity non-planar path includes three total internal 



Fig. 1. Semiconductor-diode-laser-pumped Non-Planar Ring 
Oscillator (NPRO). 


reflections and a fourth reflection at non-normal incidence from a multilayer-dielectric-coated 
curved output coupler. The total internal reflections provide reciprocal polarization rotation, and 
non-reciprocal rotation is produced by Faraday rotation when the laser material is placed in a 
magnetic field. The result is four nearly degenerate modes, two for the different eigen 
polarizations in each direction of oscillation. One of the four modes will have minimum loss. 3 
Relatively small loss differences result in single-mode unidirectional laser oscillation. Designs 
have been produced that maximize loss difference to increase the resistance to optical feedback, an 
effect that can reduce frequency stability. 5 

The beatnote lincwidth of monolithic Nd:YAG (neodymium-doped yttrium aluminum 
garnet) and Nd:GGG (neodymium-doped gadolinium gallium garnet) lasers has been observed to 
be between 3 kHz and 10 kHz for short periods of time, typically 100 ms. 6,4,7 Over longer 
periods of tim e the frequency drifts through greater excursions due to the temperature tuning rate 
which is -3.1 GHz/°C for monolithic Nd:YAG oscillators. Temperature fluctuations in the 
millidegree range produce a frequency wander of several MHz. Short term stability of a few kHz 
and MHz stability over periods of minutes, however, is adequate for a number of resonant cavity 
nonlinear frequency conversion demonstrations. And long term experiments can be carried out 
when the laser is slaved to the nonlinear cavity by active frequency control. The temperature 
sensitivity can be used for tuning, but the response is slow, typically 1 sec. Piezoelectric crystals 
bonded directly to non optical surfaces of the monolithic lasers provide limited tuning range but 
much faster response extending to frequencies of hundreds of kHz. 8 - 9 A split servo with 
feedback to both the temperature and the piezoelectric crystal can have both narrow jitter lincwidth 
and long term frequency control. 

2. Pulsed Applications 

Pulsed pumping of monolithic lasers usually results in a frequency chirp during laser output 
caused by heating. Continuously pumped lasers reach a steady-state condition with constant 
temperature. However, the continuous wave output can be used to injection seed high-power 
pulsed lasers, or the cw laser output can be time gated and amplified. Another way in which 





stable-frequency, pulsed operation is obtained from the monolithic lasers is by applying a small 
amount of modulation to the pump radiation at the relaxation oscillation frequency to produce 
spiking in the laser output. Each of these methods has been used to generate pulsed pump 
radiation for optical parametric oscillators, and the results of these investigations provide a 
heuristic review of the requirement for achieving single mode-pair oscillation of OPO's. 

a. Injection seeded Q-switched laser operation 

High-power Q-switched lasers have been a traditional source of pump radiation for singly 
resonant OPO's (SRO’s). The short buildup times and high gain of the Q-switched lasers usually 
result in highly multimode output. If the oscillation of one mode is injection seeded with a few 
mW of cw radiation, oscillation on that laser mode can buildup and deplete the population 
inversion before other laser modes have a chance to grow to significant levels. Several 
commercial lasers are available that offer injection- seeded operation with nearly time-bandwidth- 
limited output of typically 10-nsec duration. 

SRO’s have been forced to operate resonated, single-mode signal even with highly 
multimode pumping. 10 The combined frequency selection of phase matching, a dispersing 
grating, and an intracavity etalon were necessary. The multi-axial-mode pump resulted in a 
multimode output at the nonresonated idler field. When this type of frequency control is used 
with a pulse SRO pumped with an injection-seeded Q-switched laser, the result is single-mode 
output at both the signal and idler wavelengths. 11 With the addition of piezoelectric control of the 
cavity length and computer control of all the adjustable parameters, spectrographic measurements 
with 300- MHz resolution were possible. 12 These SRO’s used angle tuned LiNb0 3 pumped at 
1.06 jxm. 

The single-mode pump alone is not sufficient to produce single-mode-pair oscillation in a 
simple SRO with no frequency selection except that of phase matching. This was observed in a 
BaB 2 0 4 SRO pumped with the 354-nm third harmonic of the output of an injection seeded Q- 
switched Nd:YAG laser. 13 The single frequency pump did reduce OPO output energy fluctuation 
to 10% from 30% which was observed with multimode pumping when injection seeding of the Q- 
switched laser was blocked. The plane-parallel-cavity SRO with 1.2-cm-long BaB 2 0 4 crystal and 
3-cm overall length pumped by a 6-nsec pulse would oscillate on typically 8 axial modes at the 
signal wave. The broad tuning range of the BaB 2 0 4 OPO extending from 412 nm to 2.55 pm 
included both the 1.064- pm fundamental and 532-nm second harmonic of the laser, which were 
available for injection seeding the SRO. With seeding, the buildup time of the SRO oscillation 
was substantially decreased, and single-mode-pair oscillated was achieved. Measurements of the 
SRO threshold yielded values that were less than predicted from calculations based on reported 
values of the nonlinear coefficient of BaB 2 0 4 . This lead to a reassessment of the scale for optical 
nonlinear coefficients. 

b. Nonlinear optical coefficient measurements 

The injection-seeded spatially-filtered Q-switched Nd:YAG laser provided a suitable pump for 
phase-matched second-harmonic measurements of the nonlinear optical coefficient. High optical 


quality in both the pump beam and the nonlinear material is needed for precise control of phase 
matching which is also a requirement of these measurements. The measurements included careful 
characterization of phase matching for each observation. In addition it was determined that there 
were no extraneous spatial, temporal, or spectral components of the pump. 

The results of careful nonlinear optical coefficient measurements 14 are reproduced in 
Table I. The investigation, which at first was motivated by the BaB 2 0 4 SRO results, grew to 
include a total of six materials. This was necessary because the earlier BaB 2 0 4 measurements 
were relative to KDP, and a controversy exists over the value of the nonlinear coefficient of that 
material. The measurements described here were performed in a two-beam setup which permitted 
both relative measurements between two crystal and absolute measurements with harmonic power 
measured in one beam and fundamental power measured in the other. The reproducibility of these 
measurements was ±4% which should also be the relative accuracy. The absolute accuracy is 


estimated to be better than ±10%. 

The KDP coefficient listed in Table I 
agrees within experimental error with other 
phase-matched harmonic measurements. 15 The 
coefficients tabulated here yield a ratio 
d 22 (BaB 2 0 4 )/d3 6 (KDP) of 5.8 at variance with 
the ratio of 4.1 reported earlier. The absolute 
value for the BaB 2 0 4 coefficient, however, 
more closely agrees with the SRO threshold 
observations. A more significant variance is 
with earlier parametric fluorescence 
measurement of LiI0 3 for which the nonlinear 
coefficient was 1.73 times the value reported 
here. Relative measurements between KDP and 
LiI0 3 , however, are in agreement with earlier 
results. The difference between the 
previously reported nonlinear coefficient of 
KTP and that reported here is even larger: d t{f = 
7.3 pm/V reported earlier compared to 3.18 
pm/V reported here. It is possible that some of 
these differences are due to sample to sample 
variations of the materials. Confirmation will 
come from application of these values. The 
highly coherent laser output that is becoming 
routinely available from diode-pumped solid- 
state lasers allows greater precision in the 
determination of nonlinear optical coefficients. 


TABLE I 

NONLINEAR OPTICAL COEFFICIENTS 


Crystal 

Nonlinear optical 
coefficient^) 
(10' 12 m/V) 

KDP 

d 36 = 0.38 

KD*P 

d 36 = 0.37 

LiI0 3 

1 

H 

5%Mg0:LiNb0 3 

N 

1 

BaB 2 0 4 

‘‘at - 194 
1^22 1 - 2-2 <i) 

KTP 

d eff * 3.18 
\d l5 1 = 2.6< c > 
Id^l - 3.3 


(a) Nonlinear coefficients are given for 
1064- to 532-nm second-harmonic 
generation. 

(b) Assumes that W 3 ,l « M for BaB 2 0 4 . 

(c) Using \eUA / W 15 l = 1.25 and assuming 
^24 and d |5 have the same sign. 


c. Long-pulse-pumped SROs 

Greater frequency selection is possible with longer buildup times to parametric oscillation. Pump 
pulses of 500-nsec duration were derived from a diode-pumped NPRO by gating the cw output 
followed by multipass amplification in a flashlamp pumped laser amplifier. 16 A monolithic SRO 
was pumped by the second harmonic generated by the oscillator-amplifier laser system. Pumping 
with the 532-nm harmonic allowed noncritical phase matching in the 5%Mg0:LiNb0 3 crystal. 
The purpose of this experiment was to investigate spectral narrowing and to reduce OPO 
threshold to a level approaching that which could be achieved directly by diode-pumped lasers. 

The monolithic SRO tuned from 834 to 958 nm and 1.47 to 1.2 |im when temperature was 
adjusted between 190° and 125°C. 17 Damage limitation of the MgO:LiNb0 3 SHG crystal 
required that the 5-kW output of the laser amplifier be no longer than 500 ns. Under these 
conditions 800 W of 532-nm harmonic was generated. The ring-cavity configuration of the SRO 
allowed high efficiency; up to 60% pump depletion was observed after threshold was reached. 
About 20% of the time single-mode operation of the SRO was observed. More often, however, 
simultaneous oscillation on three axial modes was observed, and occasionally as many as 8 
modes oscillated simultaneously. Thus long-buildup time alone was not sufficient to guarantee 
single-mode oscillation. 

d. Doubly resonant optical parametric oscillation. 

A further reduction in OPO threshold is obtained with the doubly resonant oscillator (DRO) 
configuration in which the OPO is resonant at both the signal and idler wavelengths. The added 
constraint of double resonance in addition to phase matching and conservation of energy, 
however, makes stable operation of the DRO difficult. 18 The frequency stability of the NPRO 
and the mechanical stability of monolithic DRO construction are useful in overcoming this 
difficulty. A doubly resonant monolithic DRO was constructed from Mg0:LiNb0 3 with broad- 
band dielectric mirrors highly reflecting near 1.06 pm coated on the crystal. 19 A ring geometry 
was formed by using reflections from two 10-mm radius-of curvature surfaces on the ends of the 
noncritically phase-matched, 1.25-cm-long crystal and a polished flat on one side for total internal 
reflection. The DRO was pumped at 532 nm by second-harmonic pulses generated from the 
relaxation oscillations of a Nd:YAG NPRO. Pulsed operation was required because the DRO 
threshold was marginally higher than could be produced by the NPRO and harmonic generator 
when operated cw. 

A 10% modulation of the diode-laser current at 325 kHz drove the Nd:YAG NPRO into 
relaxation oscillations. The 1.06-pm fundamental pulses had 260-mW peak power and were 
efficiently converted into 400-ns, 230-mW, 532-nm pulses by externally resonant second- 
harmonic generation. The buildup of parametric oscillation occupied most of the pump pulse 
duration. Overall DRO efficiency was only 7% due to the long buildup time. After threshold was 
reached 60% pump depletion occurred. The DRO tuned between 1.02 and 1.12 pm by 
adjustment of both temperature and electric field. The most remarkable aspect of the DRO 
operation was that single-mode-pair oscillation was achieved on almost every pulse. The only 



exception was when the DRO was tuned between cluster centers and either simultaneous or 
alternating output on widely spaced modes approximately 4 nm apart were observed. The 
monolithic DRO with a pulsed single-mode pump will oscillate on a single mode pair due to the 
constraint of double cavity resonance. 

3. cw Applications 

Continuous wave operation is important for narrow-bandwidth nonlinear optical frequency 
conversion. Transient conditions can be avoided with cw operation; even the time-bandwidth 
limitations of the pulsed applications described above exceed the short term frequency stability of 
diode-pumped solid-state lasers. Resonant-cavity techniques allow high-efficiency nonlinear 
optical frequency conversion at the power levels available from frequency stabilized diode- 
pumped solid-state lasers. The use of resonant cavities external to the the lasers simplifies the 
operation of the laser, but the laser and external cavity resonances must be locked together. Stable 
and efficient second-harmonic generation has been achieved with feedback control of the cavity 
resonance, and it should be possible to stabilize optical parametric oscillation with similar 
techniques. 

a. External resonant cavity harmonic generation 

The use of resonant-cavity harmonic generation external to the laser 20 allows the problems of 
laser stabilization and harmonic generation to be treated independently while high conversion 
efficiency is maintained. Second harmonic generation with only the fundamental resonated has 
been demonstrated using diode-pumped solid-state lasers. 21 It is possible to resonate both 
fundamental and harmonic 22 but with double resonance care is required to preserve the phase 
relationship between fundamental and harmonic. 23 There are several methods that can be used to 
lock the fundamental frequency and the harmonic-generation resonant cavity. One of the simplest 
is the dither and lock technique. On resonance the reflected fundamental power is minimum. The 
optical length of the cavity can be modulated using either electro-optical or piezoelectric techniques 
or the frequency of the fundamental radiation can be modulated. An error signal is developed 
using phase sensitive detection of the reflected fundamental power. Using this technique, 
Kozlovsky, et al . 21 were able to generate 29.6 mW of cw 532-nm radiation from 56 mW of 
incident fu ndame ntal in a monolithic extemally-resonant MgO:LiNb0 3 harmonic generator. 

b. cw pumped DRO 

A second monolithic MgO:LiNb0 3 DRO with higher finesse was prepared which could be 
pumped with the available cw 532-nm radiation. 24 This OPO had an observed cw threshold of 12 
mW, and single-mode-pair oscillation was again observed. Temperature and electric-field tuning 
were used to tune the output from 1007 to 1 129 nm. Dither and lock techniques were used to 
stabilize the DRO on a single mode-pair but only with limited success. The DRO would track 
frequency tuning of the pump for about 90 MHz. Pump depletion of 78% was observed at two 
times above threshold. 


Measurements of the coherence properties 25 of this DRO demonstrated that the coherence of 
the pump was reproduced in the DRO output with little additional frequency noise. Furthermore, 
the sum of the phases of the signal and idler followed the phase of the pump. At degeneracy the 
DRO output was locked in phase with the pump laser output. It was not possible to use active 
control to stabilize the DRO for these measurements, perhaps due to shielding effects caused by 
trapped charges in the MgO:LiNb0 3 . Instead, the DRO was temperature stabilized. It could be 
operated for 20 minutes at degeneracy and for about one minute on a mode-pair off degeneracy 
without a mode hop; this is long enough for measurements of coherence and to indicate the 
potential for active stabilization. 

c. Injection locking. 

The use of a low power frequency-stabilized oscillator to stabilize the output of a higher power 
oscillator is an established technique. When injecting the low-power master laser output into the 
higher power slave laser, it is necessary to isolate the master from optical feedback, provide 
spatial mode matching, and to lock the resonances of the master and slave lasers. High-power, 
diode-pumped lasers are in the developmental stage, but injection locking has been demonstrated 
on a hybrid system with a diode-pump laser used to lock the output of an arc-lamp-pumped, high- 
power laser. 26 Ring geometry was used in the high-power laser to reduced optical feedback and 
avoided spatial hole burning effects. Locking of the cavities was achieved by the FM-sideband 
technique which provided an error signal used to drive mirrors mounted on piezoelectric 
translators and control the cavity length of the slave oscillator. It was possible to maintain 
injection locking for a slave-to-master power ratio of 400:1 with a 13-W, single-frequency output 
Analysis of phase noise showed the total additional rms phase noise of the slave over that present 
on the master oscillator was less than 0.3 radians. The measurement was limited by low- 
frequency laboratory acoustics. As the cost of semiconductor diode lasers decreases with 
advances in production, diode-pumped, high-power, solid-state lasers will become common. 27 A 
diode-pumped slave laser is expected to be much less noisy due to the more stable pumping 
provided by diode lasers, higher efficiency and therefore less heating, and the elimination of 
turbulent coolant flow. Such all-solid- state lasers are good candidates to scale the output of 
frequency stabilized lasers to higher powers. 

4. Nd:YAG Laser Frequency Stabilization 

It may be helpful to begin the discussion of frequency stabilization with the definition of terms 
and a mathematical overview. The specification of frequency noise uses quantities such as 
spectral bandwidth, linear spectral noise density, and Allen variance. It is often not possible to 
compare the laser to a more stable standard, and it becomes necessary compare two lasers with 
approximately the same level of frequency stability. If the two lasers are completely independent, 
the bandwidth of the difference frequency generated when the two laser outputs interfere at a 
photodetector, the heterodyne signal, sets an upper limit on the individual laser bandwidths. If 
the two lasers are not completely independent in such a comparison, it is possible that there will 
be some co mm on mode rejection, that is some frequency noise will be identical in the two lasers 


and not be evident in the heterodyne signal. Other measurements such as those derived from the 
signals used to stabilize a laser provide information on the lower bound of noise. Measurement 
by several different techniques reduce these uncertainties. 

The spectral bandwidth is the frequency width that would be observed using a scanning 
interferometer with sufficiently narrow transmission band. The heterodyne signal of the laser 
output and a much more stable reference displaced from the laser output by a small frequency 
difference translates the laser spectral distribution to a region in which it can be measured by 
radio-frequency spectral analysis. Mathematically the bandwidth * s width at half 

maximum of the power spectrum, 

W (f) = jE*(f)E(f), (1) 

where E(f) is defined by the Fourier transform 
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E(f) = lim -L. f E{t)e anft dt. (2) 

r->~> Jj J 

V 1 -7/2 

The normalization of Eq. (2) is useful for ergodic processes. An ergotic process is one for which 
the tim e average is equal to the ensemble average. With this normalization, the power spectrum of 
an ergodic process will approach a constant distribution as the sample time T increases. 

The time dependent electric field of a phase modulated output can be expressed as 

E(t) = E o exp {-i [2jcf o t + m ) , (3) 

where <t>(t) is a time dependent phase and f 0 is a fixed frequency. A phase noise can be 
represented by <f>(t). The instantaneous frequency shift is given by 


W) = <Kt)/2n- 

A complex spectral amplitude of frequency noise is obtained from the Fourier transform 


(4) 
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where T is the duration of the sample. The linear spectral density of frequency noise is obtained 
by avera gin g the absolute value of the spectral amplitude over many samples 

s/fl = {(S*(/)S(/)>} I/2 - 

The units oiSff) are Hz/ JTfz. One method of measuring S(f) is radio frequency analysis of a 
frequency discriminant signal, for example a voltage proportional to the instantaneous frequency 
shift of the laser output 



The Schawlow-Townes limit describes the noise due only to the random addition of 
spontaneous emission to the laser oscillation. The bandwidth of frequency noise in this limit is 


^FWHM 


2 nAv 2 h v. / P , 

Cl* 


(7) 


where Av c is the width of the laser cavity resonance, h is Planck's constant, v L is the frequency of 
the laser oscillation, and P is the output power of the laser. The linear spectral density of 
frequency noise in the Schawlow-Townes limit is constant in frequency and given by 


S 


f,o 


= Av c JThvJp . 


( 8 ) 


The condition (S f 0 ) 2 « B where B is the bandwidth of the noise applies, and it is appropriate to 
use the relationship 28 


^FWHM “ *^.0* 


(9) 


Another measure of frequency stability is the two sample variance 


<j 2 (t) 
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M-l 


2(A/-l)fTi 


LKr^ 2 - 


( 10 ) 


Here y i is the i 111 frequency measurement of M successive measurements each of duration f . A 
division of the two sample variance by the frequency of the oscillation yields the Allen variance. 29 
In practice the two-sample variance can be measured using the heterodyne signal obtained from 
two lasers slightly offset in frequency. A time and interval counter can be used to measure zero 
crossings of the heterodyne signal in successive time intervals. One data set can be manipulated 
mathematically to yield the variance as a function of individual measurement duration. 


a. Laser stabilization loop 

The control loop for laser stabilization can be modeled as four elements: the laser oscillator which 
is perturbed by a noise process Sf , a discriminator that monitors the frequency fluctuations 
with a voltage responsivity D v in units of V/Hz, a servo that amplifies the error signal with a gain 
G, and an actuator that converts the amplified error signal to changes in laser frequency with the 
coefficient K in units of Hz/V. Any stage may limit stabilization, but the discriminator is 
particularly susceptible to noise in the optical signal and technical noise such as length fluctuations 
of the reference cavity. 

Diode-pumped, solid-state, monolithic lasers usually operate at the Schawlow-Townes limit 
of frequency stability above 50 kHz, and exhibit a noise roughly proportional to 1 If below that 
frequency. A piezoelectric crystal can be bonded to the monolithic laser oscillator to serve as the 
frequency actuator. A typical response coefficient for this type of actuator is K = 1 MHz/V with a 
±20-MHz dynamic range and bandwidths flat to 200 kHz. Slower control of laser frequency with 
much larger range is achieved by temperature tuning. The servo has large gain at low 
frequencies, typically 120 dB, to control the large 1// noise and unity gain in the 50- to 100-kHz 



region. The discriminator technique that will be described here is called Pound-Drever 30 locking, 
which uses the reflection of a frequency-modulated signal from a reference cavity. Ideally the 
linear spectral density of frequency noise 5 f cl under closed-loop conditions would be reduced 
from the free running value for the laser by 
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f.cl 


*^f, laser 

I 1 +KGD I ’ 
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where K is the actuator response coefficient, S is the servo gain, and D y is the discriminator 
voltage responsivity. If stabilization is limited by noise in the frequency discriminant, the 
minimum linear spectral noise density will be 
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( 12 ) 


b. Frequency discriminant 

Resonant cavities can be used for a stabilization reference both in transmission and reflection. 
The frequency modulation technique using reflection proposed by Drever, et. al 30 has several 
advantages over locking to the side of a cavity transmission band. A frequency modulation is 
applied to the laser output before it is incident on the reference cavity. The modulation frequency 
is typically 10 MHz or greater. The FM side bands are separated from the carrier by the 
modulation frequency which is greater than the optical bandwidth of the reference cavity. If the 
carrier is centered on the resonance it will be at a minimum of reflection, and if the carrier is 
slightly detuned from resonance, the amplitude of its reflection will grow, and it will be shifted in 
phase by an amount dependent on the detuning. The FM sidebands, however, will be well 
separated from the cavity resonance and therefore be nearly totally reflected with no phase change. 
The reflected light has a resulting amplitude modulation which is dependent in magnitude and sign 
of the degree of detuning of the carrier. The modulation frequency of the AM signal is the same 
as that of the FM. 


In the simplest analysis an FM signal can be approximated by a carrier and two sidebands 


sin{cu o f + a + P sin(o) m r + y)} 

* sin(a> o r +a) + ^ sin[(<u+a> m )r +a +P ] + sin[(n) o -tu m )r +a ~P +n] . (13) 

2 

Note that the sum of the phases of the two sidebands differs by n from twice the phase of the 
carrier. A similar expansion illustrates that the sum of the phases of AM sidebands is exactly 
twice the phase of the carrier. More accurately the FM signal is represented by a series with 
Bessel function coefficients 



11 


exp{i(<a o r +/3 sin coj )} = 


U 0 (/3) + '£j k (P){e 


i k 0) t 


+ (-D « 


v -i k a) t 
K m 


}]■ 


i (0 t 


(14) 


k = 1 


The complex reflection coefficient of the resonator is 
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where and r 2 are the field reflection coefficients of the first and second mirrors respectively and 
the cavity round trip phase shift is given by 0(4/) = 2tt Af / 4/fsr w l* ere 4 /fsr is the free 
spectral range of the reference cavity. Eq. (15) assumes there are no losses in the cavity. If r l = 
r 2 = r, and A f is small compared to 4 /fsr> then 
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where 7 is the finesse of the resonator. If the reflectivity given by (16) is assigned to the carrier 
of (13) and the two side bands are assigned unity reflectivity, the total reflected intensity is 
obtained by taking the square of the absolute value of the amplitude. Only the term that is 
intensity modulated at angular frequency co m is kept, and this result is normalized by dividing by 
the average value of the incident intensity to yield the normalized error signal 

1 ( 0 ) ) 

- - WF 6f / V fsR ) sin OJ . (17) 


This result is in agreement with a more detailed derivation 31 that treats a single Fourier component 
of frequency noise rather than a fixed frequency offset and uses Eqs. (14) and (15) in the 
expansion. The detailed derivation uses an incident field that is modulated at two frequencies 

E = E. exp{i ( (0 Q t + asin +/3sin (O m t )} , (18) 


with the conditions that co 0 » co m » 0)^ and a « 1. The normalized signal is 
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where Xw = — — , G — 2 7! and tan (p = —G tan (x^lG). The added detail of Eq. 
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gives information about the frequency response of the frequency discriminant. Bandwidth, 
frequency dependent phase shift, and optimum value for the depth of modulation /3 = 1.08 are 
specified by this equation. 


The FM stabilization technique offers several advantages over the transmission method of 
stabilization. The error signal for the FM technique is at a frequency well above the frequency of 
excess noise amplitude fluctuations of the laser over the shot noise limit, thus preventing these 
amplitude fluctuations from being interpreted as frequency noise. Also the FM technique does not 
have loop delays due to resonance buildup time, and it has a slightly greater sensitivity to 
frequency change than the transmission method. 

c. Laser frequency noise 

It is important to characterize the spectral density of frequency noise of the laser for the purpose of 
determining design parameters of the control loop. The frequency discriminant is useful for this 
measurement. Calibration of the frequency discriminant was performed by directing the laser 
output through an acousto-optic modulator in addition to the electro-optical modulator that was 
introducing the 10-MHz FM sidebands. The frequency of the acousto-optic modulator was in 
turn modulated with a low frequency sinusoidal input. The fixed frequency offset of the deflected 
beam had no consequence, but the varying component served as a calibrated frequency noise 
signal. The total modulation of the laser output was as described in Eq. (18). The measurement 
of laser spectral noise density was performed under closed loop conditions, and the laser noise 
was deconvolved with knowledge of the characteristics of the loop components. This 
measurement was performed using a 40-mW commercial Nd:YAG NPRO. The results 31 are 
shown in Fig. 2. Similar results have been obtained with various techniques such as Fourier 
analysis of time domain measurements, fiber delay measurements, 32 and cavity transmission 
measurements. 33 

The spectral noise density of the Nd:YAG 
NPRO has a strong \tf dependence decreasing to 
a constant value at approximately 50 kHz. The 
white noise above 50 kHz is approximately four 
times the calculated Schawlow-Townes limit of 
0.05 Hz /yHz . The observed noise density at 1 
kHz is 20 Hz /Jh z. Corresponding values at 
this frequency are typically 10 3 Hz /^Hz for dye 
lasers and nearly 10 4 Hz/^Hz for Argon ion 
lasers. The noise of the NPRO appears largely 
due to fluctuations in the output of the diode 
pump laser used for pumping. This is 
substantiated by numerical modeling which 
shows that white noise in the diode-laser output 
would transfer to the spectral noise density with 
the same frequency dependence observed. If the 
diode laser noise were reduced to the shot noise 
limi t, the numerical calculation indicates that the 
spectral noise density of the NPRO would 



Fig. 2. Measured NPRO linear spectral density 
of frequency noise (data points) compared to the 
modeled frequency dependence for pump with a 
white noise component (solid line). 




decrease to 20 Hz Hz at DC. The low noise and potential for further noise reduction make the 
NPRO an excellent laser for frequency stabilization. 

d. Frequency stabilization of the NPRO 

The FM- stabilization or Pound-Drever technique was used to stabilize two identical NPRO lasers. 
The lasers were locked to adjacent fringes of a Fabry-Perot cavity with a finesse ofF= 22,000 
and a free spectral range of A/f SR = 6.327 GHz. The frequency separation resulting from the use 
of adjacent fringes allowed independent stabilization of the two lasers relative to the same 
reference cavity. In this measurement^ ^ the two lasers were constructed from Nd:GGG 
which has a higher Verdet constant than Nd: YAG. The higher Verdet constant provided greater 
loss difference of the lasers modes and therefore higher resistance to optical feedback. 

The experimental setup for the stabilization of the two lasers is shown in Fig. 3. The FM 
frequencies imposed by the two electro-optical modulators EOl and E02 were 10.9 and 20.3 
MHz respectively. The two laser beams were combined at the beamsplitter (BS). One of the 
beams from the beamsplitter was incident on a fast photodiode (Dl) used to observe the beating of 
the laser outputs at 6.327 GHz. This signal was mixed down to 20 kHz with a precision RF 
oscillator and analyzed with an audio spectrum analyzer. The second combined beam from the 
beamsplitter was transmitted through a polarizing beamsplitter (PBS), a quarter wave plate (X/4), 
and a mode-matching lens (LI), then reflected from the reference cavity back through the lens and 
waveplate to the polarizing beam splitter where it was reflected to a second photodiode (D2). The 
signal from the photodiode was mixed with the 10.9- and 20.3-MHz signals from oscillators f j 
and f 2 . The signals from the mixers were sent through respective low pass filters and amplified in 
the servos, and finally applied to the actuators piezoelectric crystals on the lasers. The actuator 
response was K = 450 kHz/V with a bandwidth of 500 kHz and a 20-MHz dynamic range. The 
servos had 125-dB gain at DC and unity gain at 100 kHz. The discriminant response was 0.4 
V/MHz. 



Fig. 3. Schematic diagram of dual-laser locking system. 
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Fig. 4. Radio-frequency spectral analysis of heterodyne beat note 
between two stabilized NPRO's. 


The beatnote bandwidth observed with the spectrum analyzer was 2.9 Hz (Fig. 4). 
Sampling times of 1.6 sec were used. Two-sample-variance measurements were also performed 
on the heterodyne signal from the two lasers (Fig. 5). The two-sample variance decreases with 
increased sample duration as typical of white noise until a minimum of 4.3 Hz was reached for x 
= 0.2 sec. The increase after the minimum is typical of long-term drift An analysis of the shot- 
noise limit of the discriminant signal shows that each laser has a theoretical noise bandwidth of 
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Fig. 5. Two-sample root variance of heterodyne beat note between 
two stabilized NPRO’s. 




2.2 Hz. All three of these values, 2.9-Hz beatnote, 4.3-Hz minimum two-sample variance, and 
calculated 2.2-Hz shot-noise-limited laser bandwidth are in reasonable agreement 

In a subsequent experiment 31 the heterodyne beatnote linewidth was reduced to 0.33 Hz. 
The main factor allowing this was an increase in laser output power. The 2-mW Nd:GGG lasers 
were replace with 40-mW Nd:YAG NPRO lasers. The increased power and increased sensitivity 
of the Nd:YAG lasers to feedback made additional isolation necessary. An acousto-optic 
modulator was used to shift the frequency of the lasers before reflection from the reference cavity. 
This in combination with the polarizing beam splitter and quarter wave plate provided the required 
degree of optical isolation. The increased power lowered the theoretical bandwidth due to shot 
noise to 21 mHz. However, another limit to stabilization had been reached. This limit appears to 
be mechanical or acoustic noise in the reference cavity. The excess noise of the reference cavity 
will change the free spectral range as well as the frequency of the resonances. 

The work described above is on relative stabilization of the laser frequency. The absolute 
stabilization is only as good as that of the reference cavity. Work is continuing to isolate the 
reference interferometer from laboratory noise and to achieve stabilization with narrower reference 
cavity resonances. The necessity for stabilization of the reference cavity is well known. 36 
Absolute stabilization of Nd:YAG lasers is possible by locking to atomic or molecular 
absorptions. 37 He-Ne lasers have been stabilized to 50-mHz levels. 36 Solid-state lasers have 
many properties important for frequency stable operation and may eventually be stabilized as well 
or better than He-Ne lasers. 

5 , Conclusion 

Diode-pumped monolithic solid-state lasers have intrinsically narrow linewidths. The nonplanar 
ring geometry offers advantages of high-power, unidirectional, single-mode operation and greater 
resistance to optical feedback. The laser frequency can be tuned either slowly by temperature or 
rapidly with contacted piezoelectric transducers. These lasers have been used in applications of 
injection seeding high-peak-power pulsed lasers and injection locking high-power cw lasers. The 
improved coherence is important in nonlinear-optical frequency conversion. Stable externally 
resonant second-harmonic generation has been achieved using these lasers for the generation of 
pump radiation. The low-noise free-running NPRO lasers has been useful in pumping optical 
parametric oscillators allowing developments that were difficult or impossible in the past. The 
improved coherence is important in more traditional nonlinear optical applications and has resulted 
in a reevaluation of some nonlinear optical coefficients. 

The noise properties of the diode-pumped NPRO lasers are suited to frequency stabilization 
applications. Relative stabilization to a reference cavity of 0.3 Hz has been demonstrated. There 
is now a challenge to reproduce this level of frequency stability in the output of an OPO at 
arbitrarily selected frequencies. Diode-pumped, monolithic, solid-state lasers have lower free 
running noise than other lasers and potential for higher levels of stabilization. 
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The frequency- tuning and control properties of monolithic doubly resonant optical parametric oscillators are 
analyzed for stable single- mode pump radiation. Single-axial-mode operation is observed on the idler and the 
signal for both pulsed and continuous pumping. Projections are made for tuning-parameter tolerances that are 
required for maintenance of stable single-frequency oscillation. Continuous frequency tuning is possible 
through the simultaneous adjustment of two or three parameters; thus the synthesis of specific frequencies 
within the broad tuning range of the doubly resonant optical parametric oscillator is permitted. 


1. INTRODUCTION 

An analysis of the frequency-tuning properties of doubly 
resonant optical parametric oscillators (DRO’s), based on 
both experimental observations and theoretical modeling, 
is presented. Specific details in this presentation of fre- 
quency control and synthesis apply to monolithic DRO’s 
constructed from LiNb0 3 . Where possible, however, re- 
sults are given with more general applicability. The pur- 
pose is achievement of a quantitative understanding of 
the conditions required for stable single-axial-mode para- 
metric oscillation and the resulting frequency stability 
of the DRO output. Approaches to frequency synthesis 
and continuous frequency tuning that are based on the 
simultaneous adjustment of two or three tuning variables 
are described. 

The potential of optical parametric oscillators (OPO’s) 
for the generation of tunable coherent radiation was rec- 
ognized more than twenty- five years ago. 1 The complex 
tuning properties of DRO’s were also revealed in early 
demonstrations and analyses. ^ Optical parametric os- 
cillation has been discussed in detail in a number of re- 
views, w and it is a subject treated in more general terms 
in a number of books that discuss nonlinear optics. 1 Im- 
provements in the quality of nonlinear-optical materials 
and in the coherence of pump sources led to a number of 
advances in the performance of OPO’s. Using recent ex- 
perimental results obtained with stable single-mode pump 
sources and monolithic DRO’s constructed from high- 
quality LiNbOa nonlinear-optical material, we are able 
to apply and to extend the earlier analyses. 

Resonance of both the signal and the idler frequencies, 
double resonance, offers the advantage of a lower thresh- 
old for parametric oscillation than in single resonance. 
Double resonance also provides additional frequency selec- 
tivity in OPO operation. These desirable properties of 
double resonance, however, come with a considerable in- 
crease in the complexity of tuning and with more restric- 
tive tolerances on pump stability and cavity stability. 
Diode-pumped solid-state lasers provide the required 
pump-frequency stability, and monolithic cavities provide 


the required mechanical stability in the OPO. Continu- 
ous tuning is difficult in DRO’s, which typically tune 
with axial mode hops and cluster jumps over hundreds of 
axial modes. Nevertheless, using improved pump sources 
and nonlinear-optical materials coupled with multiple- 
parameter control, DRO’s can in principle be operated 
stably and tuned continuously, thus widening their range 
of aoDlications. 

DRO’s can provide highly coherent output, re- 
producing the statistical properties of the pump with 
little additional noise. This was shown theoretically by 
Graham and Haken* in a quantum-mechanical analysis of 
the DRO, and it was demonstrated in experimental mea- 
surements of the coherence properties of the DRO. The 
quantum-mechanical analysis showed that the diffusion 
of the sum of the signal and the idler wave phaseB follows 
the phase diffusion of the pump wave adiabatically. Al- 
though the phase difference of the signal and the idler 
may diffuse in an undamped manner, the statistical prop- 
erties of a DRO are basically the same as those of an ideal 
laser, A result of these properties is the addition of only a 
small amount of phase noise in the output of the DRO 
above that present in the pump. This has been confirmed 
in coherence measurements of the output of a cw DRO. 
For periods of -1 min, the free-running DRO that was not 
servo locked oscillated on a single mode pair without a 
mode hop. That the DRO did not add significant excess 
line width over that present on the pump was demon- 
strated with the measurement of beating between the 
DRO output and an independent diode-laser-pumped solid- 
state laser during the periods between mode hops. The 
beat-note linewidth was 13 kHz, which was the expected 
value for the typically 10-kHz linewidtha of the pump 
laser and the independent reference laser. Additional co- 
herence measurements showed that the signal and the 
idler were phase anticorrelated when referenced to the 
pump laser. Also, the width of the signal-idler beat note 
with the DRO near but not at degeneracy was less than 
1 kHz. The signal-idler beat note indicates the frequency 
fluctuations added to the DRO output in addition to those 
present on the pump. 
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The results of the classical stationary analysis pre- 
sented here are consistent with the earlier analyses and 
measurements. The main purpose of this presentation is 
explanation of the complex tuning properties of the DRO 
in order to permit fuller use of its remarkable coherence 
and spectral properties. The theoretical presentation of 
Section 2 begins in Subsection 2.A with a qualitative 
overview of DRO tuning. This overview is used to estab- 
lish the extensive terminology required for the discussion. 
In Subsection 2.B the threshold condition for paramet- 
ric oscillation is reviewed and recast in terms that are 
more easily adapted to tuning calculations. The theoreti- 
cal basis of frequency selection is discussed in Sub- 
section 2.C. Experimental tuning data is presented in 
Section 3. The degree to which our theoretical model de- 
scribes the observed tuning justifies some confidence in 
its use for predictive calculations in Section 4, Results 
are summarized and discussed in Section 5, Finally, the 
properties of Mg0:LiNb0 3 that are required for modeling 
the experimental data are reviewed in Appendix A. 

2. THEORY 

A. DRO Tuning Overview 

A nonlinear-optical material pumped by intense optical 
radiation at frequency w p can provide gain at two lower 
frequencies, called the signal and the idler and related by 
the conservation-of-energy condition 

(i) p * a>, + a)i . ( 1 ) 

The parametric interaction is phase dependent, and 
proper phasing is required for energy to flow from the 
from the pump field to the signal and the idler fields. 
Phase-velocity matching ensures that the relative phases 
of the three waves do not change with propagation through 
the nonlinear material. Phase matching is described by 
the wave-vector mismatch, which for the case of collinear 
propagation can be expressed by the scalar relationship 

A k = k p - k $ - ki = (n p o) p - n $ a> $ - n<o><)/c, (2) 

where k py k 9t and k, are the respective wave-vector magni- 
tudes of the pump, the signal and the idler waves, with 
corresponding indices of refraction given by n p , n g , and 
and c is the velocity of light. Useful parametric gain 
exists in the range of signal and idler frequencies for 
which | A k | s t r/i, where l is the length of the nonlin- 
ear material. The parametric gain is maximum near 
Ak = 0. Phase matching is often achieved by controlling 
the birefringence of a nonlinear crystal through tempera- 
ture or angle of propagation. 

An OPO requires feedback at either (or both) the signal 
or the idler frequencies. If there is feedback at only one 
frequency, the device is called a singly resonant oscillator. 
Doubly resonant oscillators have feedback at both the sig- 
nal and the idler frequencies. Feedback can be provided 
by placing the nonlinear material in a cavity formed by 
two external mirrors, or in the case of monolithic OPO's, 
highly reflecting coatings can be applied directly to the 
nonlinear material. Ring-cavity configurations offer the 
advantages of reduced feedback to the pump source and 
improved OPO conversion efficiency. 11 Figure 1 illus- 
trates schematically several configurations for parametric 
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Fig. 1. (a) Schematic representation of optical parametric ampli- 
fication. Optical parametric oscillators can be formed by the ad- 
dition of mirrors that are separate from the nonlinear material, as 
shown in (b). Monolithic oscillators (c) and (d), with highly 
reflecting coatings (M's) applied directly to the nonlinear material, 
offer the advantages of low loss and rigidity that are important in 
stable, single-frequency DRO operation. Ring oscillators (d) offer 
the advantages of reduced feedback and improved conversion effi- 
ciency over standing-wave oscillators. 


oscillators. Both standing-wave and ring-cavity mono- 
lithic DRO's were used for the experimental observations 
described in this paper. The tuning properties were 
quite similar, and the same model of tuning properties 
could be used for both, because the return path length dif- 
fered little from the gain path in the ring resonators. 

Phase matching is the major factor in determining 
broad tuning properties of an OPO, although cavity reso- 
nances have the major effect on details of frequency tun- 
ing. The conditions a) p » a> t + and Ak * 0 define 
phase-matching curves. The most commonly shown OPO 
phase-matching curve is the parabolalike shape for type-I 
phase matching in a birefringent crystal, for which the 
signal and the idler waves have the same polarization and 
the pump wave has the orthogonal polarization. Fig- 
ure 2(a) shows a near-degeneracy (« f — an) section of the 
temperature-tuning curve for a LiNbO* noncritically 
phase-matched OPO. Propagation is along a crystal prin- 
cipal axis in noncritical phase matching, which reduces 
dependence on propagation direction and eliminates bire- 
fringent walk-off. 

The spectral width of the parametric gain is also deter- 
mined by phase matching. A typical spectral distribution 
for single-pass gain at fixed temperature is shown in Fig. 
2(b). Doubly resonant oscillation also entails simultaneous 
signal and idler resonance. Dispersion causes different 
cavity axial-mode frequency spacings for the two waves, 
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Fig, 2. (a) Typical OPO tuning curve near the degeneracy fre- 
quency fo * where f 9 is the pump frequency. The signal 
and the idler frequencies are shown for a LiNbOj OPO as a func- 
tion of the tuning parameter, in this case temperature. For a 
fixed value of the tuning parameter, single-pass parametric gain 
exists in bands that are centered on the phase- matching wave- 
lengths, as shown in (b). DRO’s have the added constraint that 
the signal and the idler cavity resonances must coincide in satis- 
fying the condition f 9 * f, + f, which results in output at cluster 
frequencies (c). Only two or three clusters, represented by open 
lines, are located within the gain bandwidth. Usually one clus- 
ter, represented by the longest open line, dominates. 


and the simultaneous resonance condition thus occurs only 
at intervals in frequency. The regions of simultaneous 
resonance, called cluster frequencies, are indicated in Fig. 
2(c). Early DRO's were observed to oscillate on a group or 
cluster of adjacent cavity axial modes. The wavelength of 
the cluster would at first shift continuously with tuning and 
then exhibit a discontinuous jump to another cluster of 
modes. The curves of Fig. 2 are intended to illustrate some 
general properties of the frequency tuning of DRO’s. The 
curves were calculated for the 12.5-mm-long monolithic 
MgO:LiNb0 3 oscillators, pumped at 564 THz (532 nm), 
that are discussed in this paper. 

A useful device for understanding the requirement for 
simultaneous cavity resonances is a type of diagram used 
by Giordmaine and Miller 1 ; the cavity resonances near 
the oscillating signal *n d idler frequencies are plotted as a 
function of the respective frequencies, as shown schemati- 
cally in Fig. 3. The difference between signal and idler 
axial mode spacings, Sw t and &<o it respectively, is exagger- 
ated in this figure for the purpose of illustration. One 
frequency, here the signal, increases from left to right. 
The other frequency scale, the idler, is determined by the 
first scale and the conservation-of-energy condition in 
such a way that a vertical line drawn through the diagram 
will give signal and idler frequencies that satisfy Eq. (1). 
If a signal-idler resonance pair lie on a vertical line, they 
satisfy the simultaneous resonance condition. If the tem- 
perature or the dc electric field applied to the crystal is 
changed, the position of the resonances will advance along 
the scales, one to the left and the other to the right, at 
slightly different rates because of dispersion, but the 
scales will not change. If pump frequency is ch a n ged, the 


frequencies of the cavity resonances will not change, but 
one of the frequency Beales will be displaced with respect 
to the other, and the respective resonances will move with 
that scale. 

Two types of discontinuous frequency shifts are indi- 
cated in Fig. 3. One is an axial mode hop and the other is 
a cluster jump. As a tuning variable is changed, better 
coincidence in satisfying the conservation-of-energy con- 
dition is attained on adjacent signal and idler axial modes. 
It then becomes advantageous for the oscillation frequen- 
cies to hop to the adjacent modes, to one higher in fre- 
quency and the other lower. This type of discontinuous 
frequency change is referred to as a mode hop. Other fac- 
tors such as phase matching also affect the selection of the 
oscillation frequencies. As the timing variable changes, 
phase matching also changes, and at some point it is ad- 
vantageous for the oscillator to jump to the next cluster. 
This is illustrated in the schematic tuning curve of Fig. 4. 
The signal or the idler oscillation frequency progresses 
along a cluster curve in a series of mode hops until another 
cluster curve more closely approaches phase matching. 
At that point the larger discontinuous frequency change 
of a cluster jump takes the oscillation to the next cluster 
curve. Figures 3 and 4 are only schematic, with disper- 
sion greatly exaggerated. Typically there several hun- 
dred axial modes between adjacent cluster frequencies. 



Fig. 3. Diagram 1 that shows the relationship between the DRO 
signal and idler resonance frequencies and the conservation-of- 
energy condition. Signal resonances are plotted as a function of 
signal frequency w, on an ordinary linear scale, with frequency 
increasing from left to right. The idler frequency scale is deter- 
mined by that of the signal and the relationship w, - w, + 

In the display of idler resonances, therefore, frequency increases 
from right to left. A signal-idler pair that has both resonances 
centered on a vertical line are in coincidence, satisfying w, * 
Wj + a,. In general there will be some frequency mismatch A w 
for mode pair. The frequency mismatch is the frequency 
shift that is required in order far either signal or idler resonance 
to produce coincidence. The detail on the left-hand side shows 
the frequency m i s m etch Aw for a mode pair and its components 
Aw, and Awi, which are the respective frequency displacements 
from the centers of the signal and the idler cavity resonances to 
the frequencies most favorable for parametric oscillation. Dis- 
persion is exaggerated in this schematic representation. There 
are typically hundreds of cavity axial modes between the cluster 
frequencies for which Aw » 0. 




Fig. 4. Schematic representation of a detailed portion of an ide- 
alized tuning curve for a DRO. Oscillation progresses along clus- 
ter curves in discontinuous frequency changes, called axial mode 
hops, as a tuning variable is changed. A larger discontinuous 
frequency change, a cluster jump, occurs when better phase 
matching exists on an adjacent cluster curve. 


Simultaneous resonance of signal-idler mode pairs oc- 
curs as a tuning parameter is continuously adjusted. In 
general, however, coincidence is not perfect, and oscilla- 
tion of a particular mode pair depends on the degree of 
frequency matching and phase matching. The frequency 
mismatch Ao> of a signal-idler mode pair can be defined as 
the shift in frequency that is required of either the signal 
or the idler in order to bring the two resonances into coin- 
cidence to satisfy Eq. (1). It is convenient to express the 
frequency mismatch as the sum of two components: 

A<o = Aw, + A coi . (3) 

Here Au» f is the frequency shift from the peak of the sig- 
nal resonance to the signal frequency that is most favor- 
able for oscillation for that mode pair. Correspondingly, 
A a)i is the frequency shift from the peak of the idler reso- 
nance to the idler frequency that is most favorable for os- 
cillation, as illustrated in Fig. 3. The signal component is 
measured on the signal frequency scale, and the idler com- 
ponent is measured on the idler frequency scale. The di- 
rections of these scales are opposite. One increases from 
left to right, and the other is reversed, increasing from 
right to left. Therefore Aw a and A o>, appear in opposite 
directions in Fig. 3, even though they have the same sign. 
The frequency displacements of the signal and idler from 
their respective resonance peaks are discussed in detail in 
Subsection 2.C. It is useful to consider the dependence of 
the OPO threshold on frequency mismatch and phase mis- 
match first. 

B. DRO Threshold with Imperfect Signal-Idler 
Frequency Coincidence 

Even an extremely small frequency mismatch can have 
significant effects on frequency selection and threshold of 
the DRO, particularly when cavity finesse is high. The 
threshold relationship obtained here is the same as that 
derived in the quantum-mechanical analysis by Graham 
and Haken 9 and is similar to but more detailed than the 
threshold equation given by Giordmaine and Miller.* The 
result given here is in terms of classical electromagnetic 
theory and is more easily applied to the tuning analysis 


that follows. This threshold relationship is limited to 
cavities with moderate- to-infinitesimal losses. The ef- 
fect of phase and frequency mismatch on the thresholds of 
DRO’s with arbitrary cavity losses was discussed by Falk. 12 
Falk’s results are used to estimate the conditions under 
which the threshold equation used here is appropriate. 

The threshold for oscillation is obtained by setting the 
parametric gain equal to the cavity losses. The electric 
fields of the pump, the signal, and the idler waves are ex- 
pressed in terms of complex amplitude and exponentials: 

Ej(z t t) = l / 2 [Ej(z)ex p i(kjz - wjt) + c.c.], 

where the subscript j indicates signal, idler, or pump, lc is 
the wave vector magnitude, (0 is the angular frequency, z is 
the coordinate in the direction of propagation, and t is time. 
The couped equations that describe parametric amplification 
of monochromatic plane waves traveling in the z direction 
are 7 


dE 9 

-jjj- » iK 9 EpE t * exp(iAJfez), 

(4a) 

d Ei 

- 7 — = i* t E p E t • exp(iAAz), 
az 

(4b) 

<L BL 

« iKpE 9 E t exp(-iAAz), 
dz 

(4c) 


where inks units are used and k 9 * u> 9 d 9 n/(n $ c) t = 

<*>id«fr/(n<c), and k p = <tf,cf«ir/(n,c), with n 9t n it and n p 
the respective refractive indices for the signal, the idler, 
and the pump, c the velocity of light, and d^ the effective 
nonlinear-optical coefficient. The solution used here is 
derived under the assumptions that at threshold pump de- 
pletion is insignificant and that the respective changes in 
signal and pump amplitudes, A E 9 and A E it are small com- 
pared with the amplitudes. Hence E, and E> are treated 
as constants in calculating the changes, that is, 

A E 9 - i* 9 E p Et*l sinc(AA//2) (6a) 


and 


A Ef - iKiE p E 9 *l sinc(AA//2). (5b) 

The length of the nonlinear crystal is again given by l, and 
the sine function is defined by sincx = (sinx) / x. 

For the low-loss DRO’s considered here, Eqs. (5) are ade- 
quate for modeling the parametric gain. Other solutions 
to Eqs. (4) include general monochromatic plane-wave so- 
lutions 14 that permit both pump depletion and arbitrary 
changes in E t and E t and somewhat more restrictive solu- 
tions that involve no pump depletion but have arbitrarily 
large changes in E t and A solution of the latter type 
was used 11 for threshold analysis of DRO’s with arbitrary 
strength of resonance. 

The parametric gain must compensate for both a de- 
crease in amplitude and for the phase change that is due to 
propagation in the cavity. The phasor diagram shown in 
Fig. 5 helpe to illustrate this discussion. After a round- 
trip cavity transit the signal electric-field amplitude is re- 
duced by a factor (1 - a,), and the phase is shifted by an 




Fig. 5. Phaaor diagrams schematically show amplitude losses a « 
and at and phase shifts ift 9 and & » after one round-trip cavity 
transit for the signal and the idler, respectively. At thresh- 
old the increments of electric- field amplitude added by optical 
parametric amplification, AE, and A E lt must restore the origi- 
nal fields. 


angle i fa. Similarly, the idler amplitude is reduced by 
(1 - a<), and the phase is shifted by At threshold this 
change is balanced by the increments of the electric field 
A E t and A E iy added by the parametric interaction: 


(1 - a,)exp(iifi,)E, + AE, = E. 

(6a) 

and 


(1 - Oi)exp(i^()^i + m Ei . 

(6b) 

Choosing time so that the pump amplitude is real, 
E p * |£? p |, expressing the signal and idler amplitudes as 
E 9 = |J5,|exp(i<£ # ) and Ei = |2?i|exp(i<h)» and applying the 
conditions that a„ a it fa, and </r, are all small, we can 
write Eqs. (6) as 

AE. = |E,|(a, J + exp(i</>. - iy.) 

(7a) 

and 


A£ ( = |Ej|( a i J + exp(i<fo - *y<). 

(7b) 

where y, = tan'V./a.) and y, = tan~ 1 (i//,/a,). 



Substituting Eqs. (7) into Eqs. (5) results in two equa- 
tions for the complex arguments and two equations for 
the magnitudes. The relationships for the complex 
arguments, 


<t>, + <t>i = y, + ir/2 


and 


fa + fa = yi + */2 


(8a) 

(8b) 


immediately yield y 9 = y* for the stationary solution, or 


a $ a< 

Note that Eqs. (8) are consistent with the result that the sum 
of the signal and idler phases is constant when referenced to 
the phase of the pump for stable single-mode-pair 
operation. 10 


The sum of the unpumped-cavity-round-trip phase 
shifts, 


^ a <f/ $ + , (10a) 

is useful for the purpose of comparison with the results of 
Ref. 12 and for conversion to frequency mismatch. When 
Eqs. (9) and (10a) are combined, the individual phase 
shifts can be expressed in terms of the sum by 


fa = 


a. 


a, + Qi 


(10b) 


and 




at 


a, + a* 


(10c) 


The threshold equation is obtained by taking the prod- 
uct of the two equations for the magnitudes that are ob- 
tained when Eqs. (7) are substituted into Eqs. (5), with the 
result that 


[(a. 2 + fa 2 )(a? + fa 2 )) 1 * = K. Ki E 2 l 2 sinc 2 (AAZ/2) 

= f 2 / 2 sinc 2 (Afc//2) . (11) 


The quantity T 2 is the parametric gain for perfect phase 
matching, and it is proportional to the pump intensity. 7 
With Eqs. (10b) and (10c), the threshold relationship given 
by Eq. (11) can be written in the form 


r / 2 = 


a.a, 


sinc 3 (A*//2)L (a 


1 + • 


j, + a,) 2 


( 12 ) 


Figure 6 shows the DRO threshold parameter r 2 / 2 as a 
function of the phase-shift sum as given by Eq. (12) for 
two sets of cavity losses a 9 =* a* — 0.0087 and a t = a, = 
0,0033, corresponding to finesse values 9 — 360 and 
£ - 960 of the DRO*s used in the experimental measure- 
ments. The shape of the threshold curve does not change 
significantly, but the minimum value is translated toward 
zero as the losses decrease. The width of the threshold 



V (radians) 


Fig. 6. Comparison of thresholds for DAO’s with differing cavity 
finesses. Thresholds are calculated as a function of the sum of 
the cavity-round-trip phase shifts ^ with Eq. (12) for two DRO’s 
with cavity finessa 9 9 ■ 9i • 360 and 9, ■ 9$ m 960. The 
shape of the curve does not change, but the width, defined as the 
region over which threshold is less than twice its minimum value, 
decreases for higher finesse. 



curves, defined as the region over which threshold is less 
than twice the minimum value, therefore decreases as 
cavity losses decrease. 

In the application of Eq. (12) thresholds are expressed in 
terms of cavity finesse because it is easier to measure the 
ratio of resonance width to spacing than it is to measure 
losses directly. A comparison of the above analysis with 
that of a parallel plate interferometer 14 or of optical cavi- 
ties in general 16 shows that the amplitude loss coefficients 
a, and a, are related to the cavity finesse at signal and 
idler frequencies and 9 i% respectively, by. 

SF f ** ir/a f , ■» v/cii , (13a, b) 

It is also more convenient to use frequency mismatch than 
phase shift. The components of the frequency mismatch 
are related to the phase shifts by 


Aw, 


8(o t *l/ t 
2 IT 


Aw, = 


Scojtpj 
27 r 


(14a, b) 


When Eqs. (3), (10) and (12)-(14) are combined, the thres- 
hold equation becomes 


sinc 2 (AAi/2) 


1 + 


/ 2Aw^,3 ? , \ 2 

+ • (15) 


The threshold relationship given by Eqs. (12) or (15) 
agrees with other threshold expressions under appropri- 
ate conditions. This result was obtained with a first- 
order plane-wave approximation for parametric gain. In 
the case of perfect phase matching, AA = 0, and no fre- 
quency mismatch, Aw = 0, these equations reduce to 
F 2 / 2 =* a $ a it which is the result obtained directly for this 
case. 7 Focusing and coupling to cavity modes 18 must be 
considered for quantitative threshold calculations. The 
plane-wave derivation of threshold is adequate for the 
analysis of tuning when it is necessary to know only the 
relative dependence of threshold on AA and Aw. 

There is also agreement with the central result of Falk's 
analysis 12 in the limit of high cavity finesse, that is, 
a, « 1 and a* « 1, Rewritten in the notation used 
here, Falk’s Eq. (9) becomes 



(16) 


where, defining R t = 1 - a, and Ri * 1 - O/, 

p m - a.*) „ r. i - Ri 1 

B ~ R? - R. 1 ’ R t 1 - R .* 

Equation (16) is a more accurate approximation of the 
DRO threshold applicable for arbitrary cavity loss. How- 
ever, it is unwieldy and must be evaluated as a limit when 
^ = 0 or when a, * o*. Evaluation of Eq. (12) yields 
threshold values that differ from those obtained from 
Eq. (16) by approximately the fraction (a,Oi) w * There is 
a fortuitous partial compensation for this disparity in 
the approximation for finesse given in Eqs. (13). For 
9 § ■ ■ 5, the difference between Eqs. (15) and (16) is 


less than 13% at the frequency mismatch for which 
threshold is twice its minimum value, and the difference 
decreases with decreasing frequency mismatch. For 
- 10, the difference is 4% at twice minimum 
threshold, and the agreement again improves as minimum 
threshold is approached. 

Equations (12) or (15) could be used directly to deter- 
mine the mode pair with the lowest threshold for oscil- 
lation. It is more convenient, however, to restrict the 
possible mode pairs on which oscillation may take place to 
a small number, based on frequency mismatch and wave- 
vector mismatch considerations. This is done in Subsec- 
tion 2.C, where it is shown that there are three mode pairs 
in the phase-matching bandwidth for which the frequency 
mismatch is a minimum. Which of these three mode 
pairs has the lowest threshold depends on the respective 
values of A A, Aw, and the cavity finesse. 

C. Frequency Selection in the DRO 
The selection of signal and idler frequencies in a cw DRO 
operating on a single mode pair is determined by two con- 
ditions: the conservation of energy stated in Eq. (1) and 

the minimum threshold for oscillation. An approxima- 
tion for the threshold condition was given in Eq. (15). 
The conservation-of-energy condition becomes implicit in 
the analysis of the condition of minimum threshold. In 
this analysis it is convenient to follow the approach used 
by Boyd and Ashkin 3 and to define the signal and the idler 
axial mode numbers 


to, » 2ln t /X, ■■ ln t o»,/(nc) 

(17a) 

and 


TO* * 2Zn,/A, = lnifA>i/{irc) t 

(17b) 

which are continuous variables that take on integer values 
at cavity resonances. The free-space wavelengths of the 
signal and the idler are given by A, and A,, respectively, 
and 2 ln 9 and 2 are the respective optical lengths for a 
round-trip cavity transit. The free spectral ranges or 
mode s pacings of the signal-idler resonances, 5w, and 8w it 
respectively, are the frequency changes that change the 
mode numbers to, and to, by one; that is, 4 

dm. 1 I dn,\ . , 

~ — = — In, + at,- — 1 “ Sat. 
dw, ire \ da>,f 

(18a) 

and 


dmj l( 3n,\ 

dtoj irc\ da>i/ 

(18b) 

The sum of the mode numbers 


to * to, + mi , 

(19) 


which is also a continuous variable, is useful for the de- 
scription of cluster effects. A signal frequency and an idler 
frequency that satisfy Eq. (1) and for which m is an integer 
comprise a cluster frequency pair. In general, cavity reso- 
nances are not located precisely at the cluster frequencies. 
Only the sum to, + to< must be an integer at the cluster 
frequencies; the individual mode numbers in general dif- 



fer from integers by amounts equal in magnitude but op- 
posite in sign. The cavity resonance pairs that most 
closely satisfy the conservation-of-energy condition and 
therefore that are most favorable for oscillation are also 
the resonances for which m is most yearly an integer. 
Equivalently, oscillation frequencies of a DRO are dis- 
placed from a frequency pair at which m is an integer by 
no more than one half of the respective axial mode spac- 
ings, whereas there are typically hundreds of modes be- 
tween adjacent signal or idler cluster frequencies at which 
m is an integer. 

Two further quantities that are useful for the descrip- 
tion of the mode hops and cluster jumps of DRO tuning, 
Am and A m,, are obtained by subtracting the integer 
nearest the mode number from the mode number: 


Am = m - ROUND(m) 

(20a) 

Am, = m, - ROUND(m,) . 

(20b) 


These quantities are used in the calculation of oscillation 
frequencies and tuning-variable tolerances. 

At optimum operating conditions, the quantities Ak, Am, 
and Am s will all be zero, indicating perfect phase matching 
and simultaneous cavity resonances at the desired signal 
and idler frequencies. Adjustment of three independent 
parameters is necessary in order to reach this condition. 
The discussion presented here is given in general terms 
with quantities Ak, m, and m s and in specific terms of the 
tuning or control parameters that are used in the exper- 
imental observations. The experimental observations use 
temperature T and applied potential V as adjustable 
parameters to control the output signal frequency (o r Pump 
frequency (O p is used as the required third adjustable 
parameter for the calculations. Simple Taylor's expansions 
for Ak, m, and m s were found to be adequate for modeling 
the observed frequency tuning: 


A second-order derivative is used for signal frequency be- 
cause the first-order derivatives (BAk/Ba),) mp and (3Am/ 
become zero at degeneracy, and dispersion of 
Bmjdo), is essential to the analysis. 

The notation of a partial derivative in parentheses with 
a parameter subscript to the right-hand parenthesis indi- 
cates that the parameter of the subscript is held constant 
for the differentiation. The conservation-of-energy con- 
dition is introduced through this device. Consider a 
function that is dependent on the signal, the idler, and the 
pump frequencies f * /*(w«» <*>*,<*>,), and require that the 
conservation-of-energy condition <a, + a>i = co p hold. 
Differentiation with respect to a>,, with u) p held constant, 
requires that as a>, is increased must decrease, or 

[JL\ 

da), da>i 

and, similarly, 

J[\ 

Ba) p j mt Ba) p Ba)i 

The derivatives used in Eqs. (21)— (23) are expanded in 
Table 1. The differentiation is straightforward and can 
be verified by inspection of Eqs. (2), (17), and (19). 


Table 1. Derivatives Used to Calculate Tuning 
of a Monolithic DRO 



a * - - •-> - 


dA* ,r 

+ — V + A^o, 


( 21 ) 


+ (S)j Wp ” w * #) + {T ” Tt) + 


+ m 0 , 
( 22 ) 


and 


dm,, „ , 1 d 2 ™#, \2 

m, = - — {a>, - cu^o) + 7“ — ~ 

Bo), l oo), 

+ ^ (T - To) + ^ + m ^ 



dA k Uj, toy Mi to, mi to| 

If m c*r~caT~caT 

d&k oigtog ai, to, bit to* 

1 V m cdV~cW~cW 


(23) 


L Phase-Matching Curve and Cluster Curves 
A number of equations considered below are identical ex- 
cept for the exchange of the tuning variables. An econ- 
on^y of notation is possible through the use of a general 
tuning parameter which is T, V t or a) p in the specific 
example or, more generally, is any single parameter used 
to tune a DRO. A phase-matching curve gives the signal 
frequency for which Ak - 0 as a function of the tuning 
parameter £ and as used here is denoted by tu s>PM (C). 
This curve is obtained from Eq. (21) or from a similar 
equation by setting the other adjustable parameters to 
fixed values and setting Ak = 0. The cluster curves 
give the signal cluster frequencies as a function of 
the timing parameter (. The cluster curves can be ob- 
tained from Eq. (22) or from a similar equation by setting 
m to integer values and again setting the other adjustable 
parameters to fixed values. For type-I phase matching, 
the phase-matching curve is a parabolalike curve, as 
shown in Fig. 2, and the cluster curves are a family of 
parabolalike curves. The oscillating frequencies are clos- 
est to the cluster curves near the points where the cluster 
curves intersect the phase-matching curve. 

The signal frequency separation of adjacent clusters ft # 
can be obtained from the second-order approximation 


Am(cluster) =* ±1 



Away from degeneracy the first-order term dominates, 
and Eq. (24) is approximated by 



(25) 


in agreement with Ref. 4, Phase- matching limitations 
result in a gain bandwidth with half-maximum values 
at the frequencies for which Ak « ±0.886tt// * ±ir jl 
and a corresponding signal frequency full width at half- 
maximum of 

Aa>,(Gain FWHM) - | y ) ’| • (26) 

In the specific case of the monolithic DRO, for which 
derivatives are given in Table 1, ^//(dm/dai,)^ * —(dAk/ 
d<o 9 ) mp and Aw, (Gain FWHM) * 2D 5 . Since the fre- 
quency separation of the clusters is approximately one 
half of the parametric gain bandwidth, there are two or 
three clusters within the gain bandwidth. This is true 
for any DRO in which the nonlinear crystal is the only 
dispersive component and the crystal is traversed twice in 
each round-trip cavity transit but has parametric gain in 
only one direction. 


2. Oscillation Frequencies 

The oscillation frequencies are determined by phase 
matching, the center frequencies of the signal and the idler 
cavity resonances, the frequency mismatch of the reso- 
nances, the finesse values of the resonances, and the axial 
mode a pacings. To model experimental observations, we 
calculate the frequencies of the parametric oscillation by 
the following procedure. First the signal frequency for 
phase matching <o is found for specified tuning parame- 


ters with the condition Ak = 0. Next the signal cluster 
frequency that is closest to o^pu is found with the 
condition Am - 0. If the DRO cavity has only moderate 
or low finesse and if the precise oscillating frequency and 
mode hops are not of concern, these two steps are all that 
are required. The extra resolution of frequency tuning 
can be obtained using the value of A m 9 at the cluster fre- 
quency This value is called Am^ci to indicate that it 
is calculated at the cluster frequency w,, a for the specified 
tuning conditions. 

The procedure for determining the fine details of tun- 
ing is illustrated in Fig. 7. The signal cluster frequency 
o) $i Ot obtained by the steps described above, is a reference 
point from which to start. The center of the nearest 
signal resonance is displaced from by frequency 
-Am #iC i5o) f , where 8a> 9 is the signal axial mode frequency 
separation. The center of the nearest idler resonance is 
displaced by -Amj.ci^oj, from the complementary idler 
cluster frequency m where &a>t is the idler 

mode frequency separation and <o p is the pump frequency. 
Since w f(C i is a cluster frequency with Am * 0, then 
Am it a = -Am^cu which permits the frequency mismatch 
of the signal-idler mode pair to be expressed as 


A <o = Am*,a(5<i>, — 5o>,) . (27) 

Recall that the frequency mismatch is the shift in fre- 
quency of either the signal or the idler resonance that is 
necessary in order to bring the resonance pair into coinci- 
dence to satisfy the conservation-of-energy condition. 

The displacement of the signal oscillation frequency 
w t 'OM from the center of the signal cavity resonance is 
Ao) 9t and the displacement of the idler oscillation fre- 
quency from the center of the idler cavity resonance 
is A at. From the above definitions and conditions im- 
posed on round-trip cavity phase shifts for a stationary 
solution given in Eqs. (3), (10), (13), and (14), it follows that 



Fig. 7. Signal-idler resonance diagram similar to Fig. 3, expanded 
in detail to show die relationships between quantities. The signal 
and the idler cavity resonances on which oscillation occurs are 
displaced from the respective cluster frequencies <o a and 0) l Cl for 
the general case of nonzero frequency mismatch. The DRO 
oscillating frequencies <o t Q§c and ^ divide the frequency 
mismatch Aw into the component* A <o 9 and A a> r 



Ao>, * Aftj&ij, 9i/(8<o, 


(28a) 

and 

A u>i = Aft>5ft*i ^,/(5ft*, + 5o*i? ,). (28b) 

Finally, the signal oscillation frequency is given by the 
sum of the cluster frequency plus the frequency separa- 
tion of the signal cavity resonance from the cluster fre- 
quency plus the frequency shift from the center of the 
signal cavity resonance, that is, 

= fcVci - Am^ci&o* + Aw, 


every time j changes by one. Recall that m^x is the 
value of m, on the cluater curve for which Am = 0. The 
change of the tuning parameter AfHop«p«einf that corre- 
sponds to a mode hop is a quantity that is easily measured 
experimentally. Since a mode hop corresponds to a 
change of one in m,£ i, it follows that the tuning parameter 
change that corresponds to the mode-hop spacing is 



In the evaluation of dm,£\/d( it is helpful to use the 
derivative 


= ft*«,ci “ Am.ci 


5ft* ,5ft* j(SF, + fy) 
5ft* + 5 oj,^, 


(29) 


Equations (27) and (28a) are used to obtain the second 
step of Eq. (29). 

The rate of frequency change associated with a general 
tuning parameter £ could be obtained by direct differ- 
entiation of Eq. (29) or more simply by considering the 
tuning rates of the cavity resonances. The tuning rate for 
the frequency of the signal mode is -Sco s (dm s /d£) and that 
for the frequency of the idler is -&^(dmj/d£). Tte tuning 
rates of the cavity resonances are combined to yield the rate 
of change of the frequency mismatch. 


eAw dm, dmi dm 

_ = 5<u _ + 5wj __ = ( 5<u . - 5<Ul .)— + So,,-- 


The tuning rate of the signal oscillation frequency is the 
sum of the tuning rate of the signal resonance plus the 
fraction 8<o, 9i /(8<o + 5ft*^,) of dAcu/df, that is, 


fa>«.o»c 

d( 


U*H- i9 ' + 9i) *!!h). ( 30 ) 

5ft* ( $ + 561^,1 1 af v ' l) a; / 


This tuning is limited to a small range by mode hops or 
increased cavity losses as the oscillation is pulled off the 
peaks of the cavity resonances. On a broader scale, tun- 
ing progresses along a cluster curve in a series of mode 
hops. If finesse is high, it is possible that the oscillation 
jumps back and forth between adjacent cluster curves and 
also hopB from one mode pair to the next along each of the 
cluster curves. The analysis of cluster jumps requires 
that the mode-hop structures on the two or three cluster 
curves that are closest to phase matching be compared in 
order to determine which cluster curve provides condi- 
tions most favorable for oscillation. 

The tuning rates in the regions between mode hops and 
cluster jumps, which are described by Eq. (30), are 
strongly dependent of the relative values of finesse of the 
signal and the idler cavity resonances. Some caution, 
however, is required in the use of Eq. (30). For example, 
the calculation of da> s0%c ibV and da> s 0 ^ 0 fdT for DRO's 
with nearly equal signal and idler finesse involves the small 
difference of two quantities. In such situations it is 
important that the terms on the right-hand side of Eq. (30) 
be evaluated accurately for the specified operating 
conditions. 


3. Tuning Limits and Mode Hops 

Mode hops are periodic along the cluster curves, occurring 


dft*a,ci 

~1T = “If 

5 ft*, 5 ft*, dm (32) 

5ft*, — 5ft* , 8f 

The first step simply states that the cluster frequency o), tC \ 
must change with the tuning parameter ( in such a way 
that m does not change, and the second step is accom- 
plished by using E<js. (18) and (19). It is possible to ex- 
pand the derivative dm^/d ( by first using the chain rule 
of differentiation and then Eqs. (18) and (32) to obtain 



8m t ci dm, dm, dft*«,ci dm, (33) 

d{ d{ + dal, d( d{ Satt - 8o>, H 

Another useful parameter is the maximum frequency 
shift from the cluster curve that can be achieved without 
a mode hop (ouh* - w.a). This maximum is obtained 
directly when the extreme values of Am^ C i = ±1/2 are in- 
serted into Eq. (29), yielding 


1 + ^.) 

"*** “ 8b)i9 ' + ' 


(34) 


Cavity finesse can also limit the single-parameter tun- 
ing range. It follows from Eq. (15) that threshold is 
double its minimum value when the frequency mismatch 
A a* reaches the value 

$5ft*, + y,5ft*( 

Aw “ ± 2 9,9. 

The corresponding value of Am,,ci is obtained from 
Eq. (27). On substitution into Eq. (29) the maximum dis- 
placement of signal oscillation frequency from the cluster 
frequency allowed by cavity finesse is found to be 

If finesse is large the frequency displacement allowed by 
Eq. (35) becom** significantly smaller than the frequency 
displacement required for a mode hop described by 
Eq. (34). In this case parametric oscillation on the cluster 
curve closest to phase matching ceases in a region near 
the mode hop. It is then possible for the parametric oscil- 
lation frequencies to jump to an adjacent cluster curve 
that is still within the phase-matching gain bandwidth if a 
favorable coincidence of signal and idler resonances exists 
on that cluster curve. 


Plotting the mode-hop frequency limits a in** and the 
finesse frequency limits in addition to the cluster 
curve provides additional information concerning the fine 
detail of tuning. On a broader scale, it is informative to 
display curves defining the phase- matching gain band- 
width along with the phase- matching curve. It may also 
be useful to display more than one cluster curve near the 
phase- matching curve. 

An attempt was made to keep the results of this section 
general. For application to the specific case of a mono- 
lithic DRO tuned by temperature, applied potential, and 
pump frequency, the appropriate variables and derivatives 
from Table 1 are directly substituted for the terms involv- 
ing the general tuning parameter £ Evaluation of the 
derivatives for the case of monolithic DRO’a made from 
MgO:LiNb0 3 with propagation in the x direction and 
with an electric field applied in the y direction is dis- 
cussed in Appendix A. Temperature-dependent disper- 
sion, thermal expansion, the electro-optic effect, and the 
piezoelectric effect of the nonlinear-optical material are 
used in the evaluation. Results of this evaluation for ex- 
perimental conditions described in Section 3 are given 
in Table 2. 

3. EXPERIMENTAL OBSERVATIONS 
AND MODELING 

A* Experimental Conditions 

Two monolithic DRO’s, which were described previ- 
ously, 17,1 * were used in the experimental observations. 
One DRO had lower finesse and had to be pulse-pumped in 
order to achieve the higher threshold power needed for 
parametric oscillation. The higher-finesse DRO operated 
above threshold with the available continuous pumping. 
The pump source was a diode-laser-pumped nonplanar 
ring oscillator, 19,30 constructed of neodymium-doped 
yttrium aluminum garnet (Nd:YAG) with the 1064-nm 
laser output converted to 532 nm by externally resonant 
second-harmonic generation. 31 Approximately 30 mW of 
cw pump radiation was generated The laser operated in 
a single longitudinal and a single transverse mode. For 
cw operation of the laser fundamental frequency stability 
was typically 10 kHz over short periods of time. 33,23 This 
value was doubled at the second harmonic. Higher peak 
power at similar average power was obtained by driving 
the laser into relaxation oscillations by 10% amplitude 
modulation of the diode-laser output at 320 kHz. Good 
frequency stability and high optical quality of the pump 
radiation, such as that achieved with the diode-laser- 
pumped solid-state laser, are important for obtaining sta- 
ble DRO performance. 

Both monolithic DRO’s were operated with a ring-reso- 
nator configuration. They were constructed from 5% mag- 
nesium-oxide-doped lithium niobate (MgOtLiNbC^). 24 * 25 
Each of these monolithic resonators was 12.5-mm long 
with the crystal x axis in the long direction. A ring path in 
these resonators was formed by reflections from two 
multilayer dielectric coated surfaces with 10-mm radii of 
curvature and a totally internally reflecting surface. The 
centers of curvature of the spherical surfaces were on a line 
parallel to and 180 \xm inside the flat totally internally 
reflecting surface. A drawing of the monolithic DRO's is 
shown in Fig. 8. The 532-nm pump beam was mode 


Table 2. Derivatives used to Model Tuning of 
MgO:LiNbO, Monolithic DRO’s* 
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matched for collinear propagation on the segment of the 
ring path parallel to the crystal x axis. The pump beam with 
extraordinary polarization did not follow the closed path of 
the signal and the idler waves with ordinary polarization 
because of bireflection. 

Metal coatings for electric-field tuning were applied 
to the crystal surfaces perpendicular to the y axis. The 
thickness of the crystals between the electrodes was 



2.2 mm. The finesse of both DRO’s at 1.064 fim was 
measured with the Nd:YAG laser output directly, without 
second-harmonic generation. Electric-field tuning was 
used to scan the resonators through a free spectral range, 
and transmission through the resonators gave a measure 
of resonance width relative to the mode spacing. One 
DRO had a finesse of 360, and the other had finesse of 
960. The lower-finesse device had an experimentally ob- 
served threshold for cw parametric oscillation of 35 mW, 
and the higher-finesse DRO had a threshold of 12 mW. 
The pump source could produce approximately 30 mW of 
cw radiation at 532 nm. The higher-threshold OPO was 
pumped by 532-nm second harmonic, which consisted of 
400-nsec pulses with 230-mW peak power at a 320-kHz 
repetition rate. 

The output of the DRO’s was tuned by temperature 
and electric field. Noncritical phase matching in MgO: 
LiNb0 3 was achieved for degeneracy at 107°C, and as 
temperature was increased the signal and the idler wave- 
lengths separated from the 1.064-^.m degeneracy point. 
For the tuning studies the potential applied to the crystal 
was repetitively ramped at fixed temperature. Output 
wavelength measurements were repeated at incrementally 
changed temperatures. An fj 10, 1-m grating monochro- 
mator with a 600-line/mm grating was used for wave- 
length measurement. The DRO output directed into the 
monochromator consisted of a series of pulses; these 
pulses resulted either from the pulse pumping or from the 
mode hops produced by the ramped voltage with continu- 
ous pumping. The radiation transmitted by the mono- 
chromator usually consisted of a few pulses in a narrow 
spectral band that could be correlated with the potential 
applied to the DRO electrodes. A schematic representa- 
tion of the experimental setup is shown in Fig. 9. 

B. Cluster Tuning 

The tuning of the high-finesse DRO involved spectral 
jumps back and forth between cluster curves as well as 
mode hops along the cluster curves. This behavior is il- 
lustrated in Fig. 10, where DRO output is displayed for a 
small voltage range at a constant temperature. In this 
figure output is resolved on three separate cluster curves. 
The monochromator slits were opened to provide a 5-nm 
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Fig. ft Schematic representation of the setup used for DRO tun- 
ing measurements. 
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Fig. 10. Oscillograms of cw-pumped DRO output, showing simul- 
taneous output on three cluster curves. The signal displayed 
is that produced by a photodiode placed after a monochromator 
with slits adjusted for a 5-nm bandpass. Each of the oscillo- 
grams corresponds to the same portion of the ramped voltage 
applied to the DRO. The change in applied potential is indicated. 
The oscillograms differ only in the wavelength setting of the 
monochromator, indicated for the individual traces. The output 
on the central cluster dominates and is so strong that the oscillo- 
scope trace does not return to the baseline. 


transmission width, one sufficient to resolve the individual 
clusters while transmitting a number of mode hops. The 
central cluster curve with signal wavelength near 1043 nm 
dominated. Two other cluster curves fit within the phase- 
matching gain bandwidth, and output on these curves was 
observed near 1053 and 1037 nm. Competition with the 
central cluster curve, which depletes the pump wave, is 
evident in the two cluster curves to either side. 

Three adjustable parameters were used to fit Eqs. (21) 
and (22) to the observations. A temperature-offset cor- 
rection is used to fit the calculated phase-matching curve. 
There are inaccuracies in both the absolute measurement 
of temperature and in the temperature dependence of the 
dispersion relationships that make this necessary. The 
temperature adjustment was accomplished by shifting the 
data a fraction of 1°C but could have equally well been 
done by changing the parameter AA 0 in Eq. (21). Another 
fitting parameter is required because the optical lengths 
of the DRO's are not precisely known. This fit is accom- 
plished by changing the value of in Eq. (22) and has the 
effect of adjusting the placement of the cluster curves. 
The thickness of the DRO crystals is also used as a fitting 
parameter. The electrodes do not completely cover the 
surfaces, and fringing effects are not considered. In- 
stead, it is assumed that there is a uniform electric field 
in the y direction given by E, » V/t, where V is the applied 
potential and t is an effective thickness. Adjusting the 
thickness has the effect of changing the slopes of the clus- 
ter curves and voltage-tuned phase-matching curve. It is 
interesting to note that the piezoelectric effect, in addi- 
tion to the electro-optic effect, is needed to model the ob- 
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Fig. 11. Observed and calculated tuning for the pulaed-pumped 
DRO with finesse of 360. The open vertical bar* in (a) show the 
extent of tuning observed as applied potential was ramped from 0 
to 1150 V at a constant temperature. The solid curves behind 
the vertical bars are calculated phase-matching curves for the ex- 
treme voltages. Voltage tuning for three temperatures is shown 
in (b)-(d), where the bold central curves are the calculated phase- 
matching curves, and the dashed curves indicate the limits of the 
phase-matching bandwidth. The dotted curves are calculated 
duster curves, and the filled circles are observed operating points 
of the DRO. This DRO, which has only moderate finesse, exhibits 
few jumps between duster curves as the voltage is ramped. The 
data are measurement# of the applied potential for a limited 
sampling of output frequences and do not represent individual 
mode hops. 





a 


Ui 


1060 > 
1064 5 



APPLIED POTENTIAL (V) 


I 


3 

Ui 

> 

< 


E 

c 


UI 


> 

< 

£ 


phase matching 
gain bandwidth 
duster curves 
data points 

tuning range 


Fig. 12. Observed and calculated tuning for the cw-pumped 
DRO with finesse of 960. As in Fig. 11, the open bars in (a) indi- 
cate the range of tuning as voltage was ramped, in this case be- 
tween -1150 V and 0, and the solid curves behind the vertical 
bars are the calculated phnsrun itching curves for the two ex- 
treme voltages. Voltage tuning is shown for three temperatures 
in (b Md). This DRO, which has higher finesse, exhibits a num- 
ber of frequency jumps between three duster curves as voltage 
is tuned. 






served tuning. When only the electro-optic effect is used, 
the calculated voltage-tuned cluster curves are parallel to 
the voltage-tuned phase-matching curve. ? A fourth fitting 
parameter not used here is the constant m t<0 in Eq. (23). 
Adjustment of m ti o would allow for different phase shifts 
at the mirror surfaces for the signal and the idler and the 
possibility of different cavity lengths. 3 Adjustment of m $t0 
would change the calculated position of the mode hops. 

The observed and fitted tuning curves for the DRQ's are 
shown in Figs. 11 and 12. The theory is most easily ex- 
pressed in terms of frequency. Frequency therefore is 
used as the primary ordinate scale in these graphs, and 
wavelength is included as a secondary scale for reference. 
The temperature-tuned phase-matching curves are shown 
in Figs. 11(a) and 12(a). At each temperature setting a 
range of output wavelengths is obtained by voltage tun- 
ing. In most cases the observed tuning ranges cover the 
space between and extend slightly beyond the calculated 
temperature-dependent tuning curves for the extreme 
voltages. Voltage was ramped from 0 to 1150 V, and 
when the crystal was reversed from 0 to -1150 V Fig- 
ures ll(b)-ll(d) and 12(b)-12(d) show the voltage tuning 
at selected fixed temperatures. The data are the volt- 
ages at which output was observed at selected frequen- 
cies. Calculated phase- matching curves, gain-bandwidth 
curves, and cluster curves are shown for comparison. In 
some instances the data are located in lines parallel to the 
calculated cluster curves but not on the cluster curves. 
This most likely is caused by inaccuracy in temperature 
measurement, and coincidence could be obtained by choos- 
ing a different temperature calibration for each setting. 
In practice DRO tuning may provide an accurate measure- 
ment of its temperature. The DRO sensitivity to tempera- 
ture will become more apparent below when the details of 
tuning are discussed. 

Figures 11 and 12 appear to be similar in a cursory ex- 
amination; however, one aspect of tuning's dependence on 
finesse is illustrated. The lower-finesse DRO, with the 
tuning shown in Fig. 11, usually oscillates on the single 
cluster curve nearest the phase-matc h i ng curve. Some- 
times the oscillation jumps back and forth between two 
cluster curves when they are nearly equally distant from 
phase matching. The output of the higher-finesse DRO 
with the tuning shown in Fig. 12 jumps between two or 
three cluster curves. This is in agreement with theoreti- 
cal predictions that show the tuning limit imposed on the 
9 = 960 DRO by the resonance widths; that is, the finesse 
limit of tuning is reached before the mode-hop limit of 
tuning for most conditions encountered. 

The cluster curves are also dependent on temperature. 
The data displayed in Fig. 12 are interpreted to show tem- 
perature tuning at constant voltage in Fig. 13. Here the 
data points are either interpolated from measurements of 
cluster tuning with voltages both higher and lower than 
the selected voltage or extrapolated from measurements 
of the cluster curves that nearly reach the selected volt- 
age. The calculated phase matching, gain bandwidth, 
and cluster curves are again in reasonable agreement 
with observation. 

C. Axial Mode-Hop tuning 

There is good agreement between observation and the cal- 



Fig. 13. Observed and calculated tuning for the cw-pumped DRO 
as a function of temperature. The same tuning data that was 
used in Fig. 12 is used here. A fixed voltage of -200 V was 
chosen. For the cases in which oscillation on a cluster curve was 
observed at voltages both higher and lower than this voltage, fre- 
quencies were obtained by interpolation and are represented by 
filled circles. For the cases in which cluster tuning came near 
but did not reach this voltage, frequencies were obtained by ex- 
trapolation and are represented by open circles. The dotted 
curves are portions of the calculated temperature-dependent 
cluster curves. The calculated phase-matching curve is the cen- 
tral bold curve, and the dashed curves show the approximate 
gain-bandwidth limits. 


culated voltage change required to produce a mode hop. 
Observations similar to those illustrated in Fig. 10 were 
performed under various conditions. The results are 
shown in Fig. 14, in which is displayed as a 

function of detuning from degeneracy. The calculated 
line is obtained from Eqs. (31) and (33). The tuning 
parameter is voltage, and it is necessary to substitute V 
for ( in the equations and further to substitute the appro- 
priate values from Table 2, and to evaluate the derivative 
{drn/da> t ) mp as a function of signal frequency. The approxi- 
mations given in Section 4 by relations (37) and (38) also 
work well in the evaluation of Eq. (31). 

Calculations of axial- mode-hop tuning along cluster 
curves were performed for conditions that would approxi- 
mate those used to produce Fig. 10. The same fitting 
parameters were used as in Figs. 12 and 13. An operat- 
ing temperature and center voltage were chosen to give 
three cluster curves centered on phase matching at the 
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Fig. 14. Mode-hop spacing in an applied potential as a function 
of detuning from degeneracy. The dots are data points and the 
solid line ia calculated from theory. 





observed operating frequencies. This was done by ma- 
nipulation and solution of Eqs. (21) and (22). Calculated 
phase-matching, gain-bandwidth, and cluster curves in 
this region are shown in Fig. 15(a), Calculated t uning 
along the three cluster curves with the detail of mode hops 
is shown in Figs. 15(b)- 15(d). These tuning curves were 
obtained by using Eq. (22) to calculate the cluster signal 
frequency cy^ci. Eq* (23) to calculate the signal mode num- 
ber m 9 p\ at the cluster frequency, and Eq. (29) to calculate 
the signal frequency of the oscillation. The mode-hop 
frequency limits given by Eq. (34) and the finesse fre- 
quency limits given by Eq. (35) are also shown. 

When the oscillating signal frequency differs from the 
cluster signal frequency by more than the mode-hop limit, 
it is advantageous for the oscillation to shift to another 
signal-idler resonance pair. When the oscillating fre- 
quency excursion from the cluster frequency reaches the 
finesse limit, the threshold for parametric oscillation is 
double the value it had when the oscillation frequencies 
coincided with the cluster frequencies, and the thresh- 
old increases for greater excursions of the oscillation fre- 
quency from the cluster frequency. Figure 15 illustrates 
how cluster jumps can be interspersed with the mode hops 
of a single cluster curve. For the calculation presented in 
Fig. 15, the finesse limit of tuning is reached before the 
mode-hop limit is reached on the central cluster curve 
shown in Fig. 15(c), Parametric oscillation on the central 
cluster curve usually dominates, because phase matching 
is best there. When the finesse limit of the frequency 
excursion from the cluster curve is reached, however, the 
parametric oscillation on the central cluster curve decays, 
and it is possible to have oscillation build up on an adjacent 
cluster curve before oscillation can build up on the next 
mode pair of the central cluster curve. Notice that the 
mode-hop spacing that is measured by the change in the 
tuning variable is different for the three adjacent cluster 
curves of Fig. 15(b)-15(d). Also, the relative positions of 
the finesse and mode-hop limits change with detuning 
from degeneracy. 

These comparisons show quantitative agreement be- 
tween observations and the tuning theory of Section 2. 
The tuning theory describes the cluster tuning of the DRO 
as well as the effects of cavity finesse on the cluster struc- 
ture in the occurrence of cluster jumps. Also, the theory 
is able to predict the observed spacing in the tuning vari- 
able of axial mode hope on a microscopic scale of tun- 
ing. This is done with temperature-dependent dispersion, 
thermal expansion, and electro-optical and piezoelectric 
effects. Only three fitting parameters are used: a tem- 
perature calibration, which entailed the translation of a 
temperature scale by a fraction of a degree Celsius, an 
adjustment of cavity length l of less than one wavelength, 
and the use of a effective crystal thickness that compen- 
sated for the nonuniformity of the electric field inside the 
crystal. With this agreement it is reasonable to consider 
applying the theory to analysis of the DRO for optical fre- 
quency synthesis. Specifically, the analysis addresses 
conditions that are necessary to reproduce the coherence 
of a frequency-stable pump with a small degree of tunabil- 
ity at any frequency in the tuning range of the DRO. 
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Fig. 15. Detailed display of calculated DRO tuning as a function 
of applied potential for conditions that would produce output similar 
to that shown in Fig. 10. All calculations are for a fixed 
temperature of 107.540'C. Detailed calculations of tuning for three 
cluster curves are shown in 0>Hd). Here the DRO output frequency 
is indicated by the open horizontal ban. A finesse of 960 is used.. 
The slope of the continuous portions of the detailed tuning curves 
(bHd) is dependent on the relative values of finesse at the signal and 
the idler frequencies, but in all cases the slope is much less than that 
of the cluster curves. 



4. FREQUENCY SYNTHESIS AND 
TUNING-VARIABLE TOLERANCES 

Parameter tolerances and continuous frequency tuning 
are topics that can be addressed with the theory pre- 
sented above. Knowledge of tolerances is important for 
stable DRO operation and for tuning to oscillation at 
specific frequencies. Continuous frequency tuning is of 
interest in many applications. Fixed-frequency operation 
with resolution Finer than a mode hop may be required, or 
perhaps truly continuous frequency coverage may be neces- 
sary. The DRO output frequencies lie within the widths 
of cavity resonances. The extent of continuous tuning 
depends on several factors, including frequency shifts of 
the cavity modes, the spectral range over which the con- 
servation of energy condition can be satisfied while main- 
taining oscillation within a selected mode pair, and the 
spectral range over which higher net parametric gain is not 
available on another mode pair. Multiple-parameter tun- 
ing, in which two or more parameters are synchronously 
changed, is required for continuous tuning over the full 
free spectral range of the oscillator. Single-parameter 
tuning will provide frequency coverage over small regions 
that are separated by the discrete mode hops. 

It is easiest to think of tolerances for situations in which 
only a single parameter is permitted to change. In prac- 
tice there are advantages in dealing with parameter toler- 
ances in pairs. For example, voltage and temperature 
adjustments could be used to maintain stable oscillation at 
a fixed frequency. It may not be possible to control tem- 
perature to the precision that is required if voltage is 
fixed, but the lack of adequate temperature control could 
be offset by voltage control. Feedback techniques could 
be used to adjust voltage in order to maintain stable oscil- 
lation on a signal-idler mode pair even in the presence of 
temperature fluctuations that by themselves would cause 
mode hops. The change in voltage required for stable op- 
eration could be used as an error signal that would in turn 
be used to return the temperature to the desired value. 

Simultaneous adjustment of three parameters could 
also be used to tune the output frequency of the DRO. As 
an example, consider a pump frequency that is ramped in 
some specified way. The conditions required for stable 
operation on a single signal— idler mode pair could be pro- 
vided by feedback control of the potential applied for 
electro-optic and piezoelectric tuning. The tolerance re- 
quired for phase matching would be much less stringent 
than that required for stable operation on a single signal- 
idler mode pair; adequate phase matching could be main- 
tained by temperature control, based on a functional 
relationship that is dependent on the pump frequency 
and voltage required for stable operation. With two- 
parameter tuning frequency matching could be main- 
tained, but it would not be possible to m a in ta in opt imum 
phase matching. 

A. TuningVariable Tolerances 

The parameter tolerances for stable operation are deter- 
mined by the more restrictive of two conditions. Mode 
hops are avoided by operation within a range of adjust- 
ment over which higher gain does not develop on another 
signal-idler mode pair. The range of adjustment over 
which oscillation can be maintained on a mode pair may 


be limited to a smaller value by the resonance width or 
equivalently by the DRO finesse. These tolerances are 
closely related to the mode-hop spacing and spectral limits 
of tuning that were discussed above, and they can be ob- 
tained from detailed tuning curves such as those shown in 
Fig. 15(b)— 15(d) for voltage tuning. Detailed tuning 
curves for temperature tuning and pump-frequency tun- 
ing are shown in FigB. 16 and 17, respectively. 

The conditions for the calculations displayed in Figs. 16 
and 17 are the same as those used to produce Fig. 16(c). 
These conditions are 9, = = 960 for a MgO:LiNbOs 

DRO of length l = 1.25 cm, pumped at 564 THz (532 nm) 
with signal frequency near 287 THz (1043 nm). Fitting 
parameters used in these calculations, such as an effec- 
tive thickness t = 0.277 cm and a length adjustment cor- 
responding to a change in m of 0.42, are the same as those 
required to fit the experimental data in Figs. 12 and 13. 
These characteristics are carried forward to other calcu- 
lations for the purpose of providing a specific example for 
comparison. 



Fig. 16. Calculated detailed timing as a function of temperature. 
For this calculation temperature is adjusted while other parame- 
ters are held constant at values that correspond to a point near 
the center of Fig. 15(c). Here, also, the slope of the continuous 
portions of the tuning curve are dependent on the relative values 
of signal and idler finesse, and this slope is small compared with 
the slope of the cluster curve. 



Fig. 17 . Calculated detailed tuning at a function of pump fre- 
quency. For this calculation pump frequency is adjusted while 
other parameters are held constant at values that correspond to 
a point near the center of Fig. 15(c). For pump-frequency tun- 
ing with equal signal and idler cavity finesse, the slops of the 
continuous portions of the tuning curve is -0.5. Because of the 
scale necessary to display the much greater slope of the cluster 
curve, the continuous portions of the tuning curve appear to 
be horizontal 





The detailed tuning curves of FigB. 15-17 are similar in 
many respects. The mode-hop and finesse limits of fre- 
quency displacement from the duster curve are indepen- 
dent of the tuning parameter. In each of the figures the 
slope of the continuous portion of tuning between mode 
hops is much smaller than the slope of the cluster curve. 
The slopes of the continuous portions of tuning are de- 
pendent on the relative finesse of the signal and the idler 
resonances. The case of equal finesse is shown in the cal- 
culated tuning curves. The illustrated curves show that 
voltage must be held within a tolerance of approximately 
1 V, temperature within 0.0006*C, and pump frequency 
within 7 MHz for stable operation in this example. 

Analytic approximations for the parameter tolerances 
for stable operation are not limited to a specific example. 
The range over which a parameter can be changed with- 
out causing a mode hop is obtained from the condition 
Am,# ** ±1/2. Recall that is the value m, for a point 
on the associated cluster curve, and the duster curve is 
a line that gives the values of signal frequency <o 9 and a 
tuning parameter ( for which Am — 0, with the other pa- 
rameters held at fixed values. It follows that the tuning- 
parameter tolerance iB 

(36) 

The derivative in the above equation can be evaluated 
with Eq. (33). An approximation of that equation for the 
case of type-I phase matching is given by 

8*o_4m t (37) 

af So>, - 5m, af 

where Scoo is the mode spacing at degeneracy (fi>, = a, = 
(°o - ©p/2). The difference between the idler and the 
signal mode spacing can also be expanded about 
degeneracy by using Eqs. (18) to obtain 

d*m 

8<t)i — 8o)$ m 25aio*“ «(<**» “ wo)- (38) 

d<»} 9 

Equation (36) and relations (37) and (38) can be combined 
to give the desired approximation for the mode-hop pa- 
rameter tolerance, namely, 

Ajjfapt<*nu» - (39) 

The adjustable par ymotar tolerance related to cavity fi- 
nesse can be obtained from the cavity-round-trip phase- 
shift sum for which the pumping threshold is twice its 
minimum value. From Eqs. (12) and (13) this phase-shift 
sum is 

*"■-“(5 4)- (40) 

Since ^ = 2irAm, the parameter tolerance determined by 
cavity fineBse is 

4 y*q /9t + 1/9,) 

~ dm/dC 


Table 3. Calculated Single-Parameter 
Continuous Tuning Rates and Parameter 
Tolerances for Stable Operation of 
the Finesse = 960 DRO Pumped 
at 563.6 THz with Signal Frequency 287.44 THz 


Tuning 

Parameter 

C 

df i.0- 
H 

Aflioptfiknae* 

AfontokruiM 

V (Voltage) 

72 kHz/V 

±0.98 V 

±0.80 V 

T (Temperature) 

-124 MHz/*C 

±0.00064*C 

±0.00051*C 

fp (Pump 
Frequency) 

0.499 Hz/Hz 

±6.7 MHz 

±5.4 MHz 


The parameter tolerance related to mode hops, given by 
relation (39), is zero at degeneracy and increases linearly 
with detuning from degeneracy. In practice, however, op- 
eration precisely at degeneracy was stable for tens of min- 
utes with no adjustments to the DRO. 10 The parameter 
tolerance related to cavity finesse, given by relation (41), 
remains approximately constant independent of detun- 
ing from degeneracy as long as finesse remains constant. 
Calculated tuning-variable tolerances are given in Table 3 
for the MgO:LiNbO s DRO for the conditions used to gen- 
erate Figs. 15(c), 16, and 17. 

B . Single-Parameter Tuning 

The cavity resonances associated with the signal wave and 
the idler wave can have significantly different finesse. The 
oscillating frequencies of the DRO will more closely align 
with the higher-finesse cavity resonance rather than the 
complementary resonance with lower finesse and greater 
width. If the frequencies of the cavity resonances change, 
the oscillation will follow the higher-finesse resonance 
more closely, to the extent possible without a mode hop or 
cluster jump. If the pump frequency changes, the fre- 
quency of the wave oscillating on the higher-finesse cavity 
resonance will remain more nearly constant than that of the 
wave oscillating on the lower-finesse resonance. In this 
sense the higher-finesse resonance pulls the frequencies of 
oscillation more strongly. 

It has been noted by Smith 4 that the continuous tuning of 
a DRO is relatively insensitive to tuning-parameter changes 
that change the optical length of the DRO resonator, but the 
signal frequency and the idler frequency both display 
approximately one half of the change that occurs in pump 
frequency. Single parameter continuous tuning is described 
by Eq. (30), which can be rewritten with m = m s + as 

faUom 8m 9 8*i 3m 9 3mA 

d{ ” 8*,9. + ” * H /' 

In the mathematical development used here, <o 9 and cj p are 
used as independent variables, with determined by 
Eq. (1). Choosing pump frequency as the variable 
parameter { requires the substitutions 

( am A dmt 

dUpJm, 3<*i 


Afn ntofermnc* 


(41) 


and 


ami 



8m, / 3m, \ ^ Q 

d( \da> p /*. 

The signal and the idler mode spacings, 5o>, and will 
differ by only a small amount, and if 9, and 9i are nearly 
equal, the tuning rate is do>^om/d<o p m 1/2. More gener- 
ally, for differing values of finesse, the tuning rate is in 
the range 0 < da>, t09 t/da) p < 1. Even though approxi- 
mately half the pump frequency tuning will be reflected 
in signal tuning, only a relatively small spectral range will 
be covered before a mode hop is encountered. Calculated 
single-parameter tuning rates for the special case of 9, = 
9 iy corresponding to the DRO’s described in Section 3, 
are given in Table 3. The partial derivatives needed in 
this calculation were evaluated for a>.,o*c * 287.44 THz 
(1043 nm) and T 0 = 107.51°C instead of being taken from 
Table 2. As explained in subsection 2.C.2, this procedure 
is required for the evaluation of 3/ J( osc and df s q s c 

which involve the small difference in two quantities, but 
has little effect on the other values in the table. 

C. Multiple-Parameter Tuning 

It is possible to extend the continuous tuning range by 
synchronously adjusting two or three parameters. Ad- 
justing two parameters simultaneously permits the condi- 
tions Am — 0 and Am, =* 0 to be maintained, but Ak will 
change. Adjusting three parameters simultaneously per- 
mits tuned parametric oscillation while Am, Am,, and Ak 
all remain equal to zero for a specified mode pair, and 
tuning is limited only by the extent that the parameters 
can be changed. 

Generalized tuning parameters with adjustable ft and ft 
and fixed ft are used for the discussion of two-parameter 
tuning. For the specific case treated here, any permuta- 
tion of voltage, temperature, and pump frequency can be 
used for these three parameters. The conditions Am * 0 
and Am, = 0 determine a relationship between ft and ft, 
and use of this relationship permits w, and AA to be ex- 
pressed as functions of ft. These relationships can be ob- 
tained by first differentiating Eqs. (21)-{23) with respect 
to ft, yielding 

dm _ /am\ dw, dm dm dft 

dft \d(M>J mp dft 8ft 8ft dft 

dm, _ q _ dm, doi, dm, + dm, dft ^ 

dft d<a, dft 8ft 8ft dft 

and 

dA* /8A*\ da*, j 8AA ( dAk dft 

dft \ 861,/*, dft 8ft 8ft dft 


(42a) 

(42b) 

(42c) 


Only a small spectral region is being considered, and it is 
unnecessary to consider second partial derivatives with 
respect to signal frequency. The first partial derivatives, 
however, must be evaluated for the operating conditions 
that are being considered A specific case for which this 
consideration is important is frequency tuning for voltage 
and temperature adjustment, which again involves a small 
difference of terms. Equations (42a) and (42b) are solved 
for do*, /dft and dft/dft. These values are substituted 
into Eq. (42c) to yield a value for dAA/dft. Results are 
calculated for three sets of tuning parameters and are 
listed in Table 4. The conditions used for these calcula- 
tions are the same as those used for Figs. 15-17 and 
Table 3; the partial derivatives in Eqs. (42) again were 
evaluated for © j0#c = 287.44 THz (1043 nm) and T 0 = 
1 07.56 1*C. 

Two of the examples given in Table 4 are briefly noted. 
The voltage-temperature tuning mentioned earlier is of 
interest for operation at a fixed frequency. The rate of 
change of the output frequency with applied potential 
when the voltage and the temperature are changed simul- 
taneously in order to maintain Am = 0 and Am, * 0 is 
calculated to be -7.9 kHz/V. The magnitude of this 
tuning rate is about 10 times smaller than the single- 
parameter voltage tuning rate given in Table 3 and sig- 
nificantly smaller thmn the 3.5-MHz/V tuning rate of the 
cavity resonance. The reduced sensitivity is important 
in Btable-frequency operation of the DRO. Simultaneous 
pump-frequency and voltage tuning is useful for continu- 
ous coverage of the spectral region between the mode hops 
of single-parameter tuning. Calculated tuning curves for 
this case are shown in Fig. 18. The same conditions used 
in Figs. 15-17 apply again to Fig. 18. Tuning is taken to 
the limits of A k = ±ir/i in the figure. 



Fig. 18. Calculated tuning for varying the voltage and the pump 
frequency simultaneously so as to maintain Am * 0 and Am, = 0. 
The dashed line is the cluster curve of Fig. 17. Tuning limits are 
taken at the point at which A k = ±pJl in the 1 -25-cm-long crystal. 


Table 4 Calculated Value* for Two-Parameter Tuning of a Monolithic MgO : LiNbO, DRO at 

f. * 287.44 THz or A. » 1043 nm 


Adjustable Parameter* Fixed Parameter 


At, 


Ah 

Ah 


dAA 

At, 


h - f„ h - V 
h - fp, Ct-T 
h - V. h - T 


T - 107.54*C 0.510 

V * -160 V 0.511 

f - 563.6 THi -7.9 kHi/V 


1.47 x 10*’ V/Hx 2.26 x 10*‘ (rad/m)/Hz 

(-04 x 10* u )*C/Hx -6.7 x 10*' (rad/m)/Hz 

(6.4 x 10~ 4 )*C/V 0.613 (rad/m)/V 



Two methods of three-parameter tuning are presented. 
First a method is described for achieving oscillation at a 
specified frequency while satisfying the conditions AA * 
0, Am = 0, and Am, =* 0. The description is in mathe- 
matical terms, but is analogous to what could be done 
experimentally. The first step in this method is adjust- 
ment of the temperature in order to achieve phase match- 
ing for the specified frequency. This is just a matter of 
changing the temperature to the value determined by 
Eq. (21). Next the pump frequency and the temperature 
are adjusted simultaneously in order to maintain the AA = 
0 phase-matching condition and to satisfy the condition 
Am = 0. Numerically this is done by setting A A - 0 and 
Am = 0 in Eqs. (21) and (22) and solving for T and a> p 
with V and w, held constant. Next, three parameters are 
adjusted simultaneously by solving Eqs. (21)— (23) for T, 
<i> pt and V with w, again held constant and AA, Am, and 
Am, set to zero. 

In practice this mathematical procedure would be 
analogous to setting temperature to a value calculated for 
phase matching and observing the location of the cluster 
curve nearest the phase-matched signal frequency. Next, 
the temperature and the pump frequency are adjusted 
simultaneously in order to move that cluster curve to in- 
tersect the desired frequency at phase matching. At this 
point oscillation is on the resonance nearest the specified 
frequency. Finally the temperature, the pump frequency, 
and the voltage are adjusted simultaneous in order to 
bring the cavity resonance to the desired frequency while 
maintaining coincidence of the Bignal and the idler modes 
and phase matching. 

The second method of three-parameter tuning concerns 
a situation in which oscillation is achieved with optimum 
phase matching and coincidence of the modes in satisfying 
the conservation-of-energy condition. Continuous output 
frequency tuning is possible over a limited range while op- 
timum DRO operating conditions are maintained. One 
parameter can be changed arbitrarily, but the other pa- 
rameters must be changed in a prescribed manner. The 
prescription for this change is again obtained by differen- 
tiating Eqs. (21)-(23) and this time setting all total deriva- 
tives equal to zero. The pump frequency o) p is chosen as 
the independent parameter, and the differentiation yields 

dAA / aAA\ da>, /aAA\ ^ 3AA dT f aAA dV 

dci) p \ doj , J — ^ dw p ifT d cj p 3V d (o p 

(43a) 

dm _ _ Idrn \ da>, /dm\ dm _dT + dm dV 

du) p \d<i),) m ^du) p dT dfi)j dV d to p 

(43b) 

and 

dm, q dm, do*, + dm, + dm, dT ^ dm, dV (43c) 

dwp dw, d o)p dco p dT dojp dV du) p 

Again, only the first partial derivatives, evaluated in the 
region of consideration, are required. All of the partial 
derivatives are determined by material characteristics 
and the DRO configuration, resulting in three linear equa- 
tions with three unknowns, which are solved by the usual 
methods. Continuous frequency coverage can be ob- 



-20 -10 0 10 20 
PUMP FREQUENCY SHIFT (GHz) 


Fig. 19t Calculated tuning for varying the voltage, the pump fre- 
quency, and the temperature simultaneously so as to maintain 
Am * 0, Am, ■ 0, and AA ■ 0. 


tained with an incremental series of continuous frequency 
sweeps. A calculation of tuning in this manner is shown 
in Fig. 19. The extent of tuning for the individual sec- 
tions will be limited by the range of parameter adjust- 
ment. For example, there may be a maximum voltage that 
can be applied, or the extent of pump frequency tuning 
may be limited. A limit of ±1000 V was used in Fig. 19. 

The calculations of two- and three-parameter tuning 
show that a DRO can be tuned to any frequency in its 
operating range 4 with reasonable adjustment of the tuning 
parameters. Continuous timing is possible over spectral 
ranges of approximately the extent of a free spectral 
range. Complete coverage of larger spectral regions has 
to be done by scanning a series of smaller regions. The 
control of individual parameters, particularly tempera- 
ture, requires difficult tolerances. The control problem 
can be shifted to another, more easily controlled parame- 
ter, such as voltage with multiple-parameter control of the 
DRO. The degree of correction required on the second 
parameter can then be used as an error signal for control 
of the first parameter. Fortunately the oscillating fre- 
quencies of the monolithic DRO, exclusive of mode hops, 
are relatively insensitive to voltage and temperature 
changes. If mode hope and cluster jumps are avoided, the 
frequency change of the DRO is approximately one half 
the frequency change of the pump. 

5. SUMMARY 

The theory that is used to model tuning of the DRO is 
verified at many stages. The first-order threshold ap- 
proximation agrees well with more general calculations in 
the limit of low cavity loss. The theory accurately models 
observed cluster curves for two monolithic MgO:LiNbO s 
DRO ’a. The modeling includes temperature-dependent 
dispersion, thermal expansion, and the electro-optic and 
piezoelectric effects in the nonlinear material. The ef- 
fect of DRO cavity finesse on the fine details of tuning 
gives a reasonable explanation of observed cluster jumps. 
Further substantiation of the model in the fine details of 
tuning is provided by the accurate prediction of the axial- 
mode-hop rate for tuning-parameter change. 

An understanding of DRO tuning is important for con- 
trolled stable operation. Continuous tuning rates were 
calculated for single- and multiple-parameter adjustment. 



Tolerances for stable operation were estimated. The re- 
sults of these calculations will be useful for DRO de- 
sign optimization. Multiple-parameter tuning, including 
pump-frequency adjustment, will be necessary for reach- 
ing any arbitrary frequency in the OPO operating range. 
With appropriate control the DRO will be able to produce 
stable outputs with a frequency stability as good as that 
available in the pump source. 

The DRO should find application in the generation of 
stable fixed-frequency radiation. Incremental tuning in 
controlled mode hops or cluster jumps will have applica- 
tions in spectroscopy and differential absorption lidar 
(light detection and ranging). Slow, high-resolution tun- 
ing will be possible over limited frequency ranges for spec- 
troscopic applications. 

The theory presented here could easily be extended 
to DRO configurations other than monolithic devices. 
Other degrees of freedom, such as direct length control in 
a discrete-component DRO, would provide greater versatil- 
ity in operation. Independent control of signed and idler 
cavity lengths would be useful in providing greatly ex- 
tended ranges of continuous tuning. The development 
of stable DRO operation is now possible through the com- 
bination of improved nonlinear-optical materials and 
frequency- stable laser development, such as in diode- 
pumped solid-state lasers. Optical parametric oscillators 
again appear to be on the threshold of reaching a potential 
that was first understood twenty-five years ago. 

APPENDIX A. MATERIAL PROPERTIES OF 
MgO:LiNbO a RELATED TO DRO TUNING 

A. Temperature-Dependent Dispersion 
Edwards and Lawrence 26 developed temperature-depen- 
dent dispersion equations for congruently grown LiNbOs, 
based on data reported by Nelson and Mikulyak 27 and 
Smith et al. 2B They use dispersion equations of the form 

- A ' + -h ? * B ' F ~ <A1) 

where 

F = (T - T 0 )(T + T 0 + 648), (A2) 

A is wavelength in micrometers, and T is temperature in 
degrees Celsius. Coefficients for congruent LiNbOs are 
as follows: 


Ordinary Extraordinary 


A, 

4.9048 

4.6820 

A, 

0.11775 

0.09921 

A, 

0.21802 

0.21090 

A. 

0.027153 

0.021940 

B , 

2.2314 x 10'* 

6.2716 x 10 

B* 

-2.9671 x 10'* 

-4.9143 x 10 

fl» 

2.1429 x 10'* 

2.2971 x 10 

To 

24.5 

24.5 


The material used in this work is not congruent 
LiNbOs; rather, it is 5%MgO: LiNbOs. There is little re- 
fractometric data available this material. To obtain an 
approximate set of equations for 5%MgO: LiNbOs, the ex- 
traordinary index was adjusted by changing Aie from 
4.5820 to 4.56207. This has the effect of increasing the 
calculated noncritical phase-matching temperature for 
1064-632-nm second-harmonic generation from -19.4 to 
107.04°C. The measured value for the MgO-doped mate- 
rial is 107°C. 2# This modification to the congruent dis- 
persion equations accurately reproduces the observed 
tuning curve for a singly resonant OPO that was tuned be- 
tween 0.85 and 1.48 pm by varying temperature between 
122 and ISO'C. 80 This modification also predicts para- 
metric florescence pumped at 514.5 and 488 nm when the 
crystal is tuned between 100 and 450*C (Table Al). 

B. Electro-Optic Effect 

Electro-optical, piezoelectric, and thermal expansion 
characterizations of LiNbO* are reviewed by Rauber. 31 
A somewhat more extensive tabulation of electro-optical 
coefficient measurements is given by Yariv and Yeh. 32 
Their treatment of the electro-optical effect is followed 
here. The index ellipsoid in a principle coordinate system 
is given by 



When an electric field is applied, the electro-optical effect 
is described by the modified index ellipsoid 

+ r “ £ *) x * + + + + ri * £ *) 2 * 

+ 2r u E k yz + 2r u E k xz + 2r ti E k xy » 1, (A4) 

where 


Table Al. Comparison of Meanred and Calculated Values for Parametric Fluorescence in MgOLiNbO t 


514.6- nm Pump Wavelength 488-nm Pump Wavelength 


Experimental 
Temperature Setting 
CC) 

Obe. Fluorescence 
Wavelength 
(nm) 

Calc. Temperature 
CC) 

Obe. Fluorescence 
Wavelength 
(nm) 

Calc. Temperature 

CC) 

148 

777.1 

146.70 

673.3 

143.17 

198 

741.2 

196.21 

653.1 

194.93 

248 

710.8 

245.26 

634.5 

245.18 

298 

684.6 

293.38 

617.3 

294.12 

318 

674.8 

312.52 



348 

660.8 

341.50 

601.3 

341.86 

398 

639.2 

389.31 

585.7 

390.38 

448 

619.2 . 

436.60 




r,kE k 


3 


2 ',ra,E k , 

kml 


i = 1, . . . 6 . 


For point group 3m, to which LiNbOj belongs, the fol- 
lowing relationships apply: n x = n y = n 0} n % - rc t , and 



0 

-r 22 

r is 

0 

r 22 

r ia 

0 

0 

r $ a 

0 

r$ i 

0 

r 51 

0 

0 

— rj. 

0 

0 


(A5) 


There are only four independent electro-optical coeffi- 
cients. We consider only application of an electric field 
along the y axis ( E =* E y ) t which further simplifies the 
index ellipsoid to 


(£ - r - E ')*' + fc? * r - E ’V + 6 + ir " E ’ yz ° *• 


(A6) 


The presence of the y-z cross term shows that the electro- 
optical effect results in a slight rotation of the principal 
axes. This is a small effect that accounts for less than 1% 
of the refractive-index change even at the highest fields 
that are considered here; therefore this rotation is ig- 
nored. For propagation in the x direction we have 


n y 


(—2 + 



n 0 


n*r n E, 

2 


(A7) 


and 


n s - n t . 

Yariv and Yeh” list three values for r n of LiNb0 5 , mea- 
sured with a low-frequency applied electric field for vari- 
ous optical wavelengths: 

r«(633 nm) - 6.8 x 10" u m/V, 

r u (1.15 fim) « 5.4 x 10" 11 m/V, 


where / is length, Ai is change in length, T is temperature 
in degrees Celsius and 7* - 25°C is a reference tempera- 
ture. The DRO length is measured along the crystal x 
axis. One set of coefficients, 

a u = (1.54 x 10' 6 rC“ l 


and 


0u - (5.3 x 10" 9 )°C" 2 , 

apply to expansion in the x-y plane, and a second set ap- 
plies to expansion in the z direction. The spread in the 
measurements reported by Kim and Smith suggests an ac- 
curacy of —10% in the two coefficients. 

D. Piezoelectric Effect 

The direct piezoelectric effect describes the electric polar- 
ization P that results when a stress T is applied to a mate- 
rial by the relationship 

Pi-d^Tjt, (A9) 

where dy* are the piezoelectric moduli. The converse 
piezoelectric effect describes the strain S that results 
when an electric field is applied to a piezoelectric material 
by the relationship* 7 

Sjt-d&Ei. (A10) 

Summation of the repeated indices is implied in both 
of the above equations, and the moduli dy* are the same 
in both equations. The elongation of the DRO cavity is 
given by Su. the x component of the change of a vector 
that lies in the x direction. For these measurements an 
electric field is applied in the y direction, E y , and the 
modulus dm is required for the calculation of strain. The 
symmetry of the stress tensor permits the use of a con- 
tracted subscript notation in which the modulus d 2 u is ex- 
pressed as d M . The 3m symmetry of the lithium niobate 
crystal reduces the number of independent moduli to four. 
The reduced matrix of piezoelectric moduli for the point 
group 3m is 


r»(3.39 Mm) * 3.1 x 10" 13 m/V. 

These values suggest that we use r tt * 5.5 x 10" u m/V 
near the wavelength 1.06 M m , the wavelength region at 
which our 5%MgO: LiNbO* DRO was operated. The value 
6.8 x 10" 12 m/V is from a measurement reported in 
1967” Note that this measurement was made even be- 
fore the growth of congruent LiNbOj was reported.* 4,35 
We are working with an still slightly different material 
5%MgO:LiNbO >f 24,25 and caution is required in applying 
these values. 


C. Thermal Expansion 

Thermal expansion measurements to second order in tem- 
perature for LiNbOs are reported by Kim and Smith. 35 
They express the fractional change of length with the 
quadratic function 


AJ 

l 


(A 8) 



0 0 0 d i6 ~2 dn 

df% 0 dig 0 0 

d * x das 0 0 0 


(All) 


From the relationships between the moduli, it follows that 
strain is given by 


Su “ dm ^2 ~ duE 2 m ~duE y . (A12) 

The value of the piezoelectric modulus reported by Smith 
and Welsh,* 5 d« * (2.08 x 10" u ) C/N, is used. They 
identify the LiNbOi, which they used as commercially 
grown, with a Curie point of 1165°C. Sound-propagation 
measurements were used to determine the piezoelec- 
tric moduli. 

The derivatives used in Eqs. (21H23) are expanded in 
Eqs. (18) and Table 1. Evaluation for the experimental 
conditions described in Section 3 with the material prop- 
erties described above is given in Table 2. Derivatives 
involving the electro-optic and the piezoelectric effects 


<*(7* - T r ) + fi(T - T**) 2 , 



are given with respect to the voltage applied to the elec- 
trodes on the crystal surfaces perpendicular to the y axis. 
An effective crystal thickness between the electrodes, t, is 
used, and the derivatives are given by 


BAk 3 r 2* 

W = ( “‘ n ‘ W,n ‘ } 2ci’ 


(A13) 


dm _ l 

dV irct 


. , /*22 
(n,w, + + (co,n, + <o,n t ) — 


(A14) 


and 


dm, 

~dV 




(A15) 
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note added in proof 

The tuning properties of doubly resonant oscillation in 
waveguide OPO's have been observed by Suche and 
Sohler. 39 There are many similarities between the tuning 
properties described here and those reported for the 
waveguide cavities. Piskarskas et a/. 40 have reported 
cluster effects in the output of DRO's synchronously 
pumped by cw mode-locked laser output Wong 41 has 
proposed the use of DRO's for optical frequency synthesis. 
The increased photoconductivity and the trapping of 
charges in MgO:LiNb0 3 may pose some problems in the 
use of dc electric Fields with that material. 25 
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